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SECTION  3.0 

BIAS -STRESS  EFFECTS  IN  MOS  CAPACITORS 


3.1  INTRODUCTION 

In  the  course  of  studying  the  effects  of  neutron  bombardnnent 
on  MOS  capacitors,  which  are  reported  in  Section  2.0  of  this  docu- 
ment, certain  bias-stress  effects  were  observed.  Specifically,, 
biasing  a device  into  strong  inversion  for  a period  of  time  caused 
the  storage  time,  as  determined  from  capacitance  versus  time 
measurements,  to  degrade.  A subsequent  accumulation  bias 
caused  the  storage  time  to  recover  to  near  its  pre -bias -stress 
value.  This  section  describes  our  bias -stress  findings. 

Nakagiri^^  noted  the  production  of  interface  states  by  applica- 
tion of  a high  field  of  either  polarity  to  an  SiO^  film.  Jeppson^® 
investigated  negative  bias  stress  effects  in  MOS  devices  over  the 
range  25  to  125°C  at  high  fields  and  observed  creation  of  interface 
states.  Other  related  studies  have  also  been  performed . ^ ^ 
Takino^^  reported  on  the  effect  of  a bias  stress  on  C-t  curves  for 
MOS  capacitors  on  both  n-  and  p-type  substrates.  The  results 
presented  herein  are  qualitatively  similar  to  those  of  Takino  in 
that  an  inversion  bias  is  required  to  observe  storage -time  degra- 
dation; i.e.,  a negative  bias  is  needed  for  n-type  and  a positive 
bias  for  p-type.  Certain  devices  examined  here  exhibited  orders - 
of-magnitude  more  of  a bias -stress  effect  than  that  observed  by 
Takino,  whereas  other  devices  yielded  results  similar  to  his. 
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T.  Takino,  Ph.D.  thesis,  UCLA  1976. 

^ M.  Nakagiri,  Jap.  J.  Appl.  Phys . 1^,  1610  (1974). 

K.O.  Jeppson  and  C .M.Svensson,  J.  Appl.  Phys . 48.  2004  (1977). 
A.  Goetzberger  and  H.E.  Nigh,  Proc.  IEEE  1454  (1966). 
^^B.E.  Deal,  M.  Sklar,  A.S.  Grove,  and  E.H.  Snow,  J. 
Electrochem  Soc . 114,  266  (1967). 

^^A.  Goetzberger,  A.D.  Lopez,  and  R.J.  Strain,  J.  Electro- 
chem Soc.  120,  90  (1973). 

^^I.  Kobayashi,  M,  Nakahara,  and  M.  Atsumi,  Proc,  IEEE 
249  (1973). 
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3.2  EXPERIMENTAL  PROCEDURE 


Capacitance  versus  time  measurements  were  made  on  MOS 
capacitors  using  a PAR  Model  410  C-V  System  and  were  recorded 
on  an  x-y  recorder  or  a storage  oscilloscope.  To  measure  C-t, 
a voltage  step  taking  an  n-type  device  from  -5  to  -lOV  was  applied 
(+5  to  +10V  for  p-type).  Devices  studied  were  fabricated  by 
Hughes  Aircraft  and  were  from  the  same  wafers  described  in  Sec- 
tion 2. 2. 2. 8,  i.e.,  Wet  N,  Dry  N,  Wet  P,  Dry  P,  and  KA-09. 

The  metal  electrode  thickness  was  thin  (650  I)  in  the  first  four 
cases  and  thick  (-'-I  pm)  for  the  KA-09  wafer. 

The  values  of  DC  bias  applied  to  study  degradation  and  recovery 
of  storage  time  were  ±10V  and  ±30V.  To  obtain  a long  storage  time 
initially,  an  n-type  device  was  biased  at  either  +10V  or  +30V  for 
several  hours  (-10V  or  -30V  for  p-type).  The  long -storage -time 
state  is  referred  to  here  as  "upper  saturation."  To  achieve  a 
short  storage  time  in  an  n-type  device,  either  -lOV  or  -30V  was 
then  applied  for  several  hours  (+10V  or  +30V  for  p-type).  The 
short-storage-time  state  eventually  reached  is  termed  "lower 
saturation"  here.  (This  terminology  is  not  strictly  correct  since 
an  inversion  bias  was  not  applied  long  enough  to  reach  complete 
saturation.  The  term  "lower  saturation"  is  used  here  to  indicate 
that  the  rate  of  Tg  decrease  has  slowed  significantly.)  Upon  reach- 
ing lower  saturation,  the  storage-time  recovery  process  was  then 
observed  using  an  applied  bias  of  polarity  opposite  to  that  used  to 
cause  degradation. 

3.3  RESULTS 

Figure  39  shows  storage  time  (Tg)  as  a function  of  time  under 
bias  stress  at  room  temperature  for  MOS  capacitors  from  the  Wet  N 
lot.  Storage  times  in  upper  saturation  are  in  the  range  100  to  400 
sec.  A negative  bias  stress  reduces  Tg  to  the  1 to  10  msec  range 
for  the  case  of  -30V  (0. 1 to  1 sec  for  -lOV).  The  time  required  to 
reach  "lower  saturation"  is  in  the  range  200  to  1000  min.  The  right- 
hand  portion  of  Figure  39  shows  that  complete  recovery  of  Tg  can  be 
obtained  for  a positive  bias  stress.  The  recovery  process  usually 
took  about  an  order  of  magnitude  less  time  to  take  place  than  the 
degradation  process.  Data  for  p-type  devices  are  shown  in  Figure  40. 
Note  that  recovery  of  the  Dry  P sample  was  not  complete  at  the  end 
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of  the  period  during  which  a negative  bias  stress  was  applied.  De- 
vices from  lot  KA-09  exhibited  only  a slight  bias -stress  effect  com- 
pared to  the  results  shown  in  Figure  39  and  40.  Negative  bias  stress 
for  an  extended  period  of  time  reduced  Tg  for  these  devices  by  only 
a factor  of  two  to  foixr.  The  results  for  KA-09  are  similar  in  this 
regard  to  those  obtained  by  Takino. 


Lt. 


T.  Takino,  Ph.D.  thesis,  UCLA  1976. 
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Figure  40.  Storage  time  degradation  and  re 
CO  very  for  MOS  capacitors  on  a p-type 
substrate. 
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A few  general  features  of  our  bias -stress  observations  should 
be  mentioned.  The  rate  of  degradation  of  Tg  under  inversion -bias 
stress  was  observed  to  depend  on  the  amount  of  time  a device  was 
held  in  accumulation  (upper  saturation).  The  degradation  rate  was 
relatively  slow  for  a device  previously  biased  into  accumulation 
for  a relatively  long  time.  In  addition,  the  recovery  rate  was  ob- 
served to  be  relatively  slow  for  a device  held  in  lower  saturation 
with  an  inversion  bias  for  a relative  long  time.  Capacitance  versus 
voltage  measurements  were  also  made  in  upper  saturation  and  in 
lower  saturation.  Only  a slight  shift  in  flatband  voltage  (20  to 
40  mV)  was  noted  in  most  cases  for  both  p-  and  n-type  devices. 
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For  n-type  devices,  the  shift  was  in  a positive  direction,  which  is 
consistent  with  the  movement  of  positive  ions  toward  the  A1  elec- 
trode since  a negative  bias  stress  was  applied.  For  p-type  devices, 
the  shift  was  in  a negative  direction,  which  is  consistent  with  the 
movement  of  positive  ions  toward  the  Si02-Si  interface  since  a posi- 
tive bias  stress  was  applied.  It  is  possible  that  the  observed  C-V 
shifts  are  due  to  movement  of  sodium  ions  through  the  oxide . 

Figure  41  shows  results  of  Zerbst  analysis  (described  in  Sec- 
tion 2.2.2. 3)  on  C-t  curves  for  a Wet  N device  bias -stressed  at  -30V. 
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Figure  41.  Effective  generation  lifetime  and  surface  generation 
velocity  as  a function  of  storage  time  for  an  MOS  capacitor 
subjected  to  a negative  bias  stress.  Values  of  Tgg  and  s were 
obtained  by  Zerbst  analysis  of  C-t  curves. 


Shown  is  the  effective  generation  lifetime  Tgg  and  the  surface  gener- 
ation velocity  as  a function  of  storage  time.  For  long  storage  times, 
Tgg  is  long  and  is  determined  by  generation  at  centers  in  the  deple- 
tion region  bulk  (i.e.,  Tgg  -•  Tg).  As  T^  degrades  due  to  application 
of  a negative  bias  stress,  Tgg  also  decreases.  This  is  attributable 
to  the  increasing  importance  of  carrier  generation  at  centers  at  the 
Si02-Si  interface.  The  surface  generation  velocity  is  observed  to 
increase  by  orders  of  magnitude  for  a comparable  change  in  Tg  and 
Tgg  (Figvire  41).  Thus,  the  apparent  effect  of  an  inversion  bias  stress 
is  to  temporarily  create  interface  states  which  enhance  the  rate  of 
carrier  generation,  and  thereby  decrease  Tg.  These  states  can  be 
removed,  at  least  partially  and  in  most  cases  completely,  by  an  ac- 
cumulation bias  stress.  An  alternative  explanation  is  discussed  below. 

Figure  42  shows  additional  features  of  the  degradation  and 
recovery  processes.  Illustrated  are  data  for  a Wet  N sample. 

In  42(a),  C-V  curves  are  shown  for  upper  and  lower  saturation. 

It  is  seen  that  the  inversion  capacitance  is  increased  when  a device 
is  placed  in  the  lower  saturation  state  by  an  inversion -bias  stress 
(-10V).  C-t  curves  in  upper  saturation  are  shown  in  42(b).  The 
bias  stress  applied  for  a few  hundred  seconds  to  obtain  the  first 
trace  causes  Tg  to  degrade,  so  a different  C-t  curve  is  obtained 
upon  repeating  the  measurement.  A C-t  curve  in  lower  saturation 
is  shown  in  42(c).  Note  that  the  initial  and  final  values  of  capaci- 
tance are  larger  here  than  in  42(b).  Upon  applying  an  accumulation 
bias  (+10V)  to  cause  recovery  to  occur,  we  observed  that  the  C-V 
curve  recovers  before  Tg  recovers.  This  is  shown  in  42(d)  and 
42(e).  In  42(d),  the  inversion  capacitance  has  recovered  to  near 
its  upper -saturation  value  whereas  the  C-t  curve  in  42(e)  still  yields 
a short  T^.  (Note  that  the  C-V  curve  is  flat  in  inversion  in  Fig- 
ure 42(d)  because  Tg  is  very  short.  At  the  same  sweep  rate,  the 
C-V  curve  in  upper  saturation  (42(a))  is  not  flat  because  Tg  is  long; 
i.e.,  the  sweep  rate  is  such  that  a deep -depleted  condition  is 
achieved.)  Following  this,  a process  we  term  "temporary  re- 
covery" occurs,  as  illustrated  in  Figure  42(f).  It  appears  that  Tg 
is  long  from  the  first  part  of  the  C-t  curve.  However,  Tg  abruptly 
appears  to  be  short  again  near  the  end  of  the  trace.  After  a bias 
stress  on  the  order  of  20  to  100  minutes,  "permanent"  recovery 
of  Tg  is  observed.  (Data  of  the  type  shown  in  Figure  42(f)  are  only 
observed  if  the  C-t  trace  is  begun  in  deep  depletion  instead  of 
inversion. ) 


The  recovery  process  depicted  in  Figure  42(f)  was  also  ob- 
served to  take  on  another  form,  depending  on  the  time  following 
application  of  a recovery  bias  at  which  observations  were  made. 
Figure  43  shows  data  obtained  a short  time  (minutes)  after  appli- 
cation of  +30V  recovery  bias  to  an  n-type  device  (a  device  differ- 
ent from  that  of  Figure  42).  The  first  part  of  the  C-t  curve  is  very 
similar  to  that  shown  in  Figure  42(f).  However,  a second  portion 
is  then  observed  in  which  the  final  value  of  capacitance  increases. 

A trace  of  this  type  indicates  that  Tg  has  not  yet  recovered  and  also 
that  the  capacitance  has  returned  to  its  lower  saturation  condition. 

It  should  be  noted  that  the  temporary  recovery  phenomena  of  Fig- 
ures 42(f)  and  43  arise  because  the  accumulation -bias  recovery 
process  is  perturbed  by  making  a C-t  measurement.  To  make  this 
measurement,  an  inversion  bias  is  applied  to  the  device  which 
causes  degradation  to  occur  once  again.  Thus,  in  Figure  42(f) 
is  initially  long  but  is  degraded  abruptly  during  the  C-t  measure- 
ment as  a result  of  the  inversion  bias. 

3.4  DISCUSSION 

A qualitative  explanation  for  a portion  of  the  present  bias- 
stress  findings  can  be  given  based  on  the  work  of  Takino.  How- 
ever, certain  aspects  of  the  data  are  not  yet  accounted  for,  and 
further  work  will  be  required  to  achieve  the  needed  insight.  Takino 
considered  the  effect  of  positive  ions  in  the  oxide  on  the  carrier 
generation  rate  at  interface  states.  For  an  n-type  substrate,  as 
positive  ions  approach  the  interface  the  electron  emission  rate  de- 
creases. In  the  p-type  case,  the  approach  of  positive  ions  causes 
the  emission  rate  to  increase.  These  effects  are  due  to  the  cou- 
lombic  interaction  between  the  ion  and  the  interface  state. 

The  present  C-V  data  taken  in  lower  and  upper  saturation  are 
explainable  in  terms  of  positive  ion  motion  in  the  oxide.  In  ad- 
dition, if  Takino 's  model  is  correct,  then  the  C-t  (i.e.,  Tg)  data 
can  also  be  accounted  for  qualitatively  on  the  basis  of  ion  motion. 

In  n-type  devices,  for  example,  application  of  a negative-bias 
stress  would  move  the  positive  ions  away  from  the  interface, 
thereby  increasing  the  carrier  generation  rate  and  decreasing  the 
storage  time.  A positive -bias  stress  reverses  the  process. 

^°T.  Takino,  Ph.D.  thesis,  UCLA  1976. 
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Figure  43.  Capacitance  versus  time  for  an  MOS  capacitor 
illustrating  the  degradation  of  Tg  and  the  increase  in 
inversion  capacitance  that  can  occur  during  a C-t  meas- 
urement of  a device  that  had  a recovery  (accumulation) 
bias  applied  for  a short  time  (minutes). 

The  increase  in  inversion  capacitance  observed  in  lower  satu- 
ration can  be  accounted  for  qualitatively  in  terms  of  the  lateral 
motion  of  positive  ions.  Consider  a p-type  device  as  an  example.  A 
positive -bias  stress  drives  ions  perpendicularly  through  the  oxide 
to  the  interface,  with  a resulting  decrease  in  storage  time.  The 
fringing  field  present  gives  rise  to  the  laterally  depleted  region 
shown  in  Figure  1.  In  addition,  this  fringing  field  will  move  ions 
laterally  along  the  interface  beyond  the  region  under  the  gate 
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electrode.  These  ions  can  then  deplete  an  additional  portion  of  the 
semiconductor  in  the  laterial  region,  thereby  causing  the  measured 
inversion  capacitance  to  increase. 

The  above  qualitative  model  does  not  appear  to  account  for  the 
following  experimental  observations:  (1)  the  rate  of  degradation  of 
Tg  was  about  an  order  of  magnitude  slower  than  the  rate  of  recovery: 

(2)  capacitance  recovery  occurred  before  the  onset  of  storage -time 
recovery  upon  application  of  an  accumulation  bias.  One  further 
experiment  performed  near  the  end  of  this  program  also  needs  to 
be  explained.  The  temperature  of  an  n-type  (Wet  N)  MOS  capaci- 
tor with  a positive  bias  (30V)  applied  was  raised  to  either  50  or 
100°C.  Next,  the  polarity  of  the  bias  was  reversed  (-30V)  and  bias- 
temperature  stress  then  continued  for  typically  one  hour.  Upon  re- 
turning the  device  to  room  temperature  while  maintaining  the  nega- 
tive bias,  it  was  subsequently  found  that  the  time  required  to  degrade 
Tg  significantly  using  a negative -bias  stress  had  increased  by  a 
factor  of  2 to  5 compared  to  that  required  in  an  identical  device  that 
had  not  been  bias -temperature  stressed  (e.g.,  data  of  Figure  39). 

(This  experiment  was  replicated  using  several  Wet  N devices.)  In 
addition,  the  negative-bias  stress  at  elevated  temperature  did  not 
significantly  degrade  Tg  at  room  temperature;  i.e.,  the  value  of  Tg 
was  still  near  upper  saturation.  This  recent  finding  appears  to  be 
inconsistent  with  the  positive -ion  model  described  above.  However, 
it  does  appear  to  be  consistent  with  our  observation,  noted  above, 

that  the  rate  of  degradation  under  inversion -bias  stress  depends  on  * 

the  amoiint  of  time  a device  is  previously  held  in  accumulation  at 
room  temperature.  Holding  a device  in  accumulation  at  elevated 
temperature  slows  the  degradation  rate  much  more  significantly. 

The  bias -stress  observations  reported  here  are  interesting 
from  several  viewpoints,  including:  (1)  the  potential  impact  of  these 
effects  on  the  study  of  MOS  device  properties:  (2)  the  reliability  of 
MOS  devices;  (3)  the  potential  relationship  between  the  observed 
phenomena  and  interface -state  studies  at  various  laboratories  (in- 
cluding studies  of  radiation-induced  interface  states).  Additional  i 

study  of  bias -stress  effects  are  needed  before  a comprehensive 

model  can  be  developed.  Such  a model  will  be  required  before  the  . | 

importance  of  the  observed  phenomena  can  be  realistically  ! 

assessed.  i 


SECTION  4.  0 


LEAKAGE  CURRENT  PHENOMENA  IN  IRRADIATED 

SOS  DEVICES* 

4.  1 INTRODUCTION 

Silicon -on-sapphire  MOS  devices  are  being  developed  for 
radiation-hardened  electronic  systems  because  they  have  several 
advantages  over  their  MOS -on -bulk -silicon  counterparts:  reduced 
transient  photoresponse;  elimination  of  latchup;  increased  packing 
density  and  speed.  A disadvantage  of  SOS  devices  is  the  occurrence 
of  a significant  ionizing -radiation-induced  back-channel  leakage 
current  in  n-channel  transistors.  35-37  This  phenomenon  has  been 
shown  to  be  due  to  the  trapping  of  positive  charge  in  the  sapphire 
near  the  silicon-sapphire  interface  which  induces  a conducting 
channel  in  the  silicon.  \ye  have  performed  a detailed  study 

of  ionizing  radiation  effects  on  SOS  devices  with  emphasis  placed 
on  determining  the  mechanisms  of  back -channel  leakage  current 
phenomena . 

Several  related  topics  are  treated  in  this  section,  including: 

(a)  comparison  of  radiation-induced  leakage  current  for  devices 
with  wet  and  dry  gate  oxides;  (b)  radiation-induced  reduction  of 
back-channel  leakage  current;  (c)  bombarding  electron  energy  de- 
pendence of  leakage  current  production  and  leakage  current  reduc- 
tion. This  last  topic  was  explored  using  6-  to  45-keV  electrons. 

The  reader's  attention  is  particularly  directed  to  the  discussion  of 
radiation-induced  reduction  of  leakage  current.  This  novel  mech- 
anism constitutes  a potential  remedy  for  leakage  current  problems 
in  SOS  devices  employed  in  a space  radiation  environment. 

' Published  in  the  IEEE  Transactions  on  Nuclear  Science  24, 
2119  (1977). 

^^D.  Neamen,  W.  Shedd,  and  B.  Buchanan,  IEEE  Trans. 

Nucl.  Sci.  21_,  211  (Dec.  1974). 

36  R. A.  Kjar  and  J.  Peel,  IEEE  Trans.  Nucl.  Sci.  208 
(Dec.  1974). 

37k. M.  Schlesier,  IEEE  Trans.  Nucl.  Sci.  21,  152 
(Dec.  1974). 

Harari,  Appl.  Phys.  Letters^,  25  (1976). 

^E.  Harari  and  D.J.  McGreivy,  IEEE  Trans.  Electron 
Devices  1277  (1977). 
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4.2  EXPERIMENTAL  CONSIDERATIONS 


Irradiations  were  performed  at  room  temperature  using  the 
Co^*^  source  at  Northrop  Research  and  Technology  Center  and  a 
Jeol  JSM-2  scanning  electron  microscope.  The  Co^®  dose  rate 
was  4.3  rads(Si)/s  in  most  cases  but  in  some  experiments  was  in- 
creased to  64  rads(Si)/s.  Silicon-on-sapphire  devices  examined 
in  this  study  were  obtained  from  several  sources.  Most  of  the  in- 
formation presented  herein  resulted  from  studies  on  custom  n-channel 
SOS  transistors  fabricated  at  Hughes  Aircraft  (HAC).  Sapphire 
wafers  containing  a silicon  epitaxial  layer  (initial  thickness  0.5  pm) 
were  obtained  by  Hughes  from  Union  Carbide  and  devices  containing 
both  wet  and  dry  gate  oxides  were  then  fabricated.  All  processing 
steps  were  identical  for  these  devices  except  the  gate  oxidation. 

Wet  oxides  ('^980  A --  pyrogenic)  were  grown  at  925°C  whereas  the 
dry-oxide  (~930  A)  growth  temperature  was  1000°C.  The  gate  met- 
allization on  most  of  the  transistors  was  purposely  kept  thin  (750 

e 

±100  A)  to  facilitate  studies  with  a low-energy  electron  beam.  Chan- 
nel width  was  8.0  mils  and  channel  length  was  0.21  ±0.02  mil.  We 
also  investigated  radiation -induced  leakage  current  in  n-channel 
transistors  contained  in  radiation-hardened  CMOS /SOS  inverters 
obtained  from  three  manufacturers. 

Back-channel  leakage  current  I^  was  obtained  from  measure- 
ments of  drain  current  versus  gate  voltage  for  a constant  value  of 
Vjs>  typically  5V.  Such  characteristics  were  recorded  on  an  x-y 
plotter  by  applying  a ramp  voltage  to  the  gate  electrode  and  moni- 
toring drain  current  by  use  of  a current-to -voltage  converter.  In 
addition  to  obtaining  I^  from  I^js^^gs  characteristics,  we  also 
used  a circuit  to  switch  momentarily  to  measurement  bias  condi- 
tions when  it  was  necessary  to  minimize  the  time  during  which  ir- 
radiation bias  conditions  were  interrupted.  A constant  I^jg  (equal 
to  I^)  was  achieved  only  if  Vgg  was  several  volts  negative;  con- 
duction along  the  front  channel  and  at  island  edges  must  be  turned 
off  before  back-channel  leakage  current  alone  is  observed.  Thresh- 
old voltage  was  determined  for  SOS  transistors  from  measurements 
of  versus  Vgg  in  the  conventional  way.  Front-channel  field- 

effect  mobility  was  also  determined  from  such  measurements.  In 
addition,  drain  characteristics  were  used  to  obtain  front-channel 
conductance  mobility. 
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4.3  re:sults 

4.3.  1 Threshold  Voltage  Shifts  and  Channel  Mobility  Degradation 

Threshold  voltage  was  measured  as  a function  of  dose  up  to 
10*^  rads(Si)  (Co^*^  --  64  rads(Si)/s)  for  dry-  and  wet-oxide  custom 
transistors.  During  irradiation  and  measurement,  V^g  was  equal 
to  5V,  and  irradiations  were  performed  with  two  values  of  Vggi  0 
(grounded  gate)  and  5V.  Both  device  types  were  found  to  be  quite 
radiation  tolerant.  Under  worst-case  biasing  conditions,  the  maxi- 
mum value  of  threshold  voltage  shift  observed  up  to  10^  rads  was 
1 • 3V.  Interface -state  effects  were  evident  (i.e.,  positive  shift  in 
threshold  voltage)  at  high  doses  for  wet-oxide  transistors  but  not 
for  dry-oxide  units.  McLean  et  al.  previously  observed  such 
effects  to  be  more  important  in  wet-oxide  devices. 

Degradation  of  front- channel  mobility  was  found  to  be  negli- 
gible for  dry-oxide  transistors  after  irradiation  to  10^  rads.  How- 
ever, in  wet-oxide  devices  after  the  same  total  dose,  mobility  de- 
graded by  roughly  a factor  of  two  and  this  reduction  is  qualitatively 
consistent  with  the  interface -state  effects  noted  for  such  units  in 
threshold  voltage  measurements.  That  is,  increased  scattering  is 
expected  when  a significant  number  of  interface  states  have  been 
introduced  by  radiation. 

4.3.2  Leakage  Currents  in  Irradiated  Wet-  and  Dry-Oxide  Devices 

Figure  44  shows  leakage  current  versus  dose  for  wet-  and  dry- 
oxide  transistors  irradiated  using  Co^®  at  a dose  rate  of  4.3  rads 
(Si)/s.  (Di  scussion  of  the  reduction  of  shown  in  the  right-hand 
side  of  this  figure  is  given  in  the  next  subsection.  ) The  irradiation 
bias  values  were  Vgg  = 0 and  V^g  = 5V,  To  measure  I^,  Vgg  was 
switched  momentarily  to  about  -7V.  Leakage  currents  for  dry- 
oxide  devices  are  observed  to  be  substantially  larger  than  for  their 
wet  coxinterparts . Behavior  is  similar  in  both  cases  in  that  plateau 
regions  are  evident  at  high  doses  and  a significant  increase  in 
occurs  at  intermediate  doses.  (Similar  data  have  previously  been 
obtained  by  Harari . ^9) . ^ slight  decrease  in  leakage  current 

23 

F.B.  McLean,  H.E.  Boesch,  Jr.,  and  J.M.  McGarrity, 

IEEE  Trans.  Nucl.  Sci.  23.  1506  (1976). 

^®E.  Harari,  Appl.  Phys . Letters  2^,  25  (1976). 

Harari  and  D.J.  McGreivy,  IEEE  Trans.  Electron 
Devices  ZA,  1277  (1977). 
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BACK-CHANNEL  LEAKAGE  CURRENT  (A) 


DOSE  IradsISDI 

Figure  44.  Back-channel  leakage  current  vs  dose  for  dry-  and  wet- 
oxide  n-channel  SOS  transistors  irradiated  with  V„g  = 0 and 
V(jg  = 5V  at  a dose  rate  of  4.3  rads(Si)/s.  The  reduction  in  leak- 
age current  in  the  right-hand  portion  of  this  figure  was  accom- 
plished by  continuing  to  irradiate  beyond  10^  rads  with  V^jg  = 

Vgs  = (Note  the  change  in  the  dose  scale.) 

(~20%)  is  evident  as  the  dose  is  increased  from  10^  to  10^  rads.  By 
analogy  to  effects  at  the  Si02-Si  interface,  this  decrease  may  be  due 
to  the  introduction  of  negatively  charged  states  at  the  Si-Al202  inter- 
face which  would  compensate  for  a portion  of  the  positive  charge 
trapped  in  the  sapphire.  (For  data  obtained  at  a dose  rate  of  64  rads 
(Si)/s,  this  phenomenon  was  evident  but  only  at  about  the  5%  level. ) 

At  the  Northrop  Co^®  facility,  a device  to  be  irradiated  is  placed 
in  the  exposure  room  and  the  radiation  source  is  then  raised  up  into 
the  room.  This  takes  approximately  5 minutes,  during  which  time 
the  dose  rate  is  time  varying.  A maximum  rate  is  reached  when  the 
Co^®  source  is  completely  raised.  As  discussed  below,  it  was  found 
that  a small  photocurrent  dominated  behavior  shown  in  Figure  44  for 
the  wet-oxide  device  for  doses  up  to  200  rads  (the  approximate  dose 
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received  during  raising  of  the  Co  source  for  the  case  of  a maxi- 
mum rate  of  4.3  rads(Si)/s).  This  current  scaled  with  dose  rate 
and  served  to  mask  true  back-channel  behavior  at  low  doses.  To 
examine  only  the  radiation -induced  back -channel  leakage  compo- 
nent of  measured  current,  both  the  photocurrent  and  the  pre- 
irradiation  current  were  subtracted  from  measured  values. 

Figure  45  shows  radiation-induced  back-channel  leakage  current 
versus  dose  based  on  data  in  Figure  44.  The  large  vertical  error 
bar  for  the  lowest-dose  wet-oxide  data  point  reflects  primarily  the 
accuracy  of  the  photocurrent  correction  made  at  that  dose.  Also 
shown  in  Figure  45  (solid  circles)  is  a horizontal  translation  of 


10“  10'  lO’  lO’  10'  10^  io‘ 
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Figure  45.  Radiation-induced  leakage  current  vs  dose  for  wet- 

and  dry-oxide  SOS  transistors  irradiated  at  a rate  of  4.3  rads 
(Si)/s.  These  data  were  obtained  by  correcting  the  data  of 
Figure  44  for  pre  -irradiation  current  and,  in  the  wet-oxide 
case,  for  a photocurrent  component. 
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■wet-oxide  data  by  a factor  of  90  in  dose.  This  translation  illustrates 
that  both  'wet  and  dry  data  can  be  described  by  the  same  po'wer  la'w 
(current  “ (dose)^*^)  over  approximately  three  orders  of  magnitude 
variation  in  radiation-induced  leakage  current.  If  we  consider  wet 
and  dry  data  separately,  the  (dose)l*5  dependence  in  the  wet  case  is 
observed  over  a leakage  current  range  of  approximately  two  orders 
of  magnitude  for  doses  between  500  and  10^  rads  (Si).  For  dry  de- 
■vices,  this  dependence  is  also  noted  over  approximately  two  orders 
of  magnitude  variation  in  current  for  doses  between  20  and  500  rads. 

Before  interpreting  the  results  of  Figure  45,  it  is  useful  to  con- 
sider the  temperature  dependence  data  of  Figure  46.  Shown  is  leak- 
age current  versus  1000/T  for  wet-  and  dry-oxide  transistors  be- 
fore and  after  irradiation  to  several  doses.  Before  irradiation,  data 
for  wet  devices  exhibit  two  regimes  with  activation  energies  differing 
by  a factor  ot  two.  The  high -temperature  regime  is  most  likely  due 
to  thermal  generation  in  the  drain  depletion  region.  For  a single 
generation  level  at  mid-gap,  the  dominant  term  in  the  expression  for 
thermal  generation  rate  will  be  n^,  the  intrinsic  carrier  concentra- 
tion. Over  the  temperature  range  of  interest,  n^  is  approximately 
proportional  to  exp(-0.60  eV/kT)  and  thus  the  0.57  eV  slope  in  Fig- 
ure 46  is  consistent  with  the  domination  of  generation  current  at  high 
temperatures.  Harari  and  McGreivy  obtained  con'vincing  evidence 
that  pre -irradiation  leakage  current  at  room  temperature  in  wet- 
oxide  devices  very  similar  to  those  studied  here  is  attributable  to 
thermal  generation.  Although  the  present  results  are  consistent  with 
their  findings,  it  is  interesting  to  note  in  Figure  46  that  this  statement 
only  applies  at  room  temperature  and  above.  At  lower  temperatures 
for  unirradiated  wet  de-vices  and  for  both  de-vice  types  under  all  other 
conditions  shown  in  the  figure,  measured  leakage  current  is  dominated 
by  back-channel  phenomena. 

Formation  of  a channel  in  the  silicon  adjacent  to  the  Si-Al203 
interface  occurs  due  to  buildup  of  positive  trapped  charge  in  the 
sapphire.  The  quantity  induces  a carrier  density  " 

suiting  leakage  current  is  given  by  e F where  F is  applied 

field,  A is  cross-sectional  area,  and  is  back-channel  carrier 
mobility.  In  general,  and  will  be  related  nonlinearly.  Ad- 
ditionally, p^jj,  is  expected  to  vary  with  Work  of  Brews, 

Harari  and  D.J,  McGreivy,  IEEE  Trans.  Electron 
Devices  1277  (1977). 

Brews,  J.  Appl.  Phys.  46,  2181  (1975);  _46,  2193  (1975). 
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Figure  46.  Leakage  current  vs  reciprocal  temperature  for  wet- 
and  dry-oxide  SOS  transistors.  Data  obtained  before  irradi- 
ation and  after  irradiation  to  several  values  of  total  dose  are 
shown. 
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Chen  and  Muller,  and  Guzev  et  al.  has  revealed  that  front- 
channel  mobility  in  an  MOS  transistor  in  weak  inversion  is  constant 
for  very  low  values  of  n^j^^  (<10^  cm"^)  and  increases  with 
for  higher  values  of  this  quantity  (but  still  in  weak  inversion). 

Brews'  analysis  indicates  that  this  will  occur  when  front-channel 
mobility  is  dominated  at  low  carrier  concentrations  by  a relatively 
high,  inhomogeneously  distributed  fixed  interface  charge  density. 

As  the  gate  voltage  is  increased,  thereby  increasing  in  the 

front  channel,  screening  will  increase  with  a resultant  reduction 
in  carrier -density  fluctuations  . Hence,  mobility  increases . (This 
reasoning  can  also  be  applied  to  the  back  channel  in  an  SOS  device 
but  in  this  case  n^^  is  increased  by  increasing  through  ioniz- 
ing radiation.  ) Beyond  weak  inversion,  scattering  dominates  and 
mobility  decreases  with  increasing  inversion-layer  carrier 
density.  40-42 

The  temperature  dependence  of  shown  in  Figure  46  yields 
activation  energies  (E^^)  which  are  roughly  equal  except  for  the  case 
of  post-irradiation  dry-oxide  data.  (We  exclude  the  thermal  genera- 
tion regime  from  this  discussion.)  Chen  and  Muller^^  observed  that 
front-channel  mobility  in  MOS  transistors  decreased  with  decreasing 
temperature  in  the  weak  inversion  regime  and  plots  of  such  data 
yielded  an  activation  energy  of  ~0.  1 eV.  Data  in  Figure  46  are  quali- 
tatively consistent  with  their  findings  in  that  is  also  observed  to 
decrease  with  decreasing  temperature,  which  suggests  that  the 
temperature  dependence  of  leakage  current  is  attributable  to  that 
of  back-channel  mobility.  Chen  and  Muller  proposed  a temperature - 
activated  front-channel  mobility  at  low  inversion -layer  carrier  den- 
sities, with  being  an  effective  barrier  associated  with  a non- 
uniform  potential  distribution  at  the  SiO^-Si  interface.  (Analysis 
by  Brews 40  indicates  that  interpreting  channel  mobility  data  in  terms 
of  temperature  activation  may  be  oversimplified.)  This  effective 
barrier  is  lowered  by  an  increase  in  and  thus  one  should  observe 

an  that  decreases  with  increasing  in  the  weak  inversion 

regime.  The  dry-oxide  data  in  Figure  46  are  qualitatively  consistent 
with  this  concept  in  that  E^  decreased  due  to  irradiation  (i.e., 
higher  back-channel  n^j^^).  The  observation  that  wet-oxide  devices 

40 

J.R.  Brews,  J.  Appl.  Phys . 46,  2181  (1975);  46,  2193  (1975). 
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yielded  approximately  the  same  value  of  E_  under  all  conditions 

Oi 

shown  in  this  figure  is  qualitatively  consistent  with  the  front- 
channel  data  of  Chen  and  Muller,  where,  at  a given  temperature, 
mobility  (and  thus  E^)  was  observed  to  be  nearly  constant  in  very 
weak  inversion  for  several  temperatures  near  300  K. 

Pre -irradiation  data  in  Figure  46  for  the  dry-oxide  device  are 
very  similar  to  post-irradiation  wet-oxide  findings.  This  suggests 
that  a back  channel  already  exists  in  dry  devices  before  irradiation, 
whereas  in  the  wet  case  pre -irradiation  leakage  current  at  room 
temperature  is  dominated  by  thermal  generation  in  the  drain  deple- 
tion region.  Harari  and  McGreivy^^  have  explained  such  behavior 
in  terms  of  the  existence  of  a processing -produced  n-type  layer  in 
the  silicon  adjacent  to  the  Si-Al203  interface  for  dry-oxide  devices. 
They  ruled  out  trapped  charge  in  the  AI2O3  as  being  responsible 
for  pre -irradiation  leakage  current  in  such  devices  because  they 
were  unable  to  reduce  by  optical  excitation.  , 

Interpretation  of  the  room-temperature  data  of  Figure  45  is 
now  given.  The  same  functional  dependence  of  radiation-induced 
leakage  current  on  ionizing  dose  is  exhibited  for  wet-  and  dry- 
oxide  devices,  which  suggests  that  the  same  physical  processes 
are  occurring  in  both  cases.  The  factor -of-ninety  displacement 
in  dose  exhibited  in  the  figure  (i.e.,  a larger  dose  required  to 
achieve  a given  value  of  in  wet  devices)  then  needs  to  be  ac- 
counted for.  At  some  constant  value  of  radiation-induced  leakage 
current,  the  product  of  mobility  and  carrier  density  in  the  back 
channel  for  the  dry-oxide  case  must  be  equal  to  that  for  wet  oxide. 
Such  a comparison  can  be  made  for  leakage  currents  ranging  from 
~10"^  to  ~3  X 10"'^A  and  for  doses  ranging  from  ~10  to  ~200  rads 
for  dry  oxide  and  from  ~10^  to  •'-Z  x 10“^  rads  for  wet  oxide.  Since 
the  rate  at  which  the  mobility-carrier  density  product  increases 
with  dose  is  the  same  for  wet  and  dry  devices,  one  is  tempted  to 
conclude  that  at  a given  value  of  leakage  current  mobilities  are 
equal  and  carrier  densities  are  equal  for  the  two  transistor  types. 
The  displacement  in  dose  exhibited  in  Figure  45  would  then  be  at- 
tributed to  a larger  density  of  traps  in  dry-oxide  devices  than  in 
their  wet-oxide  counterparts.  The  fact  that  leakage -current  sat- 
uration occurs  at  approximately  the  same  dose  suggests  that  the 
same  trap  is  responsible  for  observed  behavior  in  both  cases. 
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Acceptance  of  a significantly  larger  trap  density  in  dry-oxide  de- 
vices then  suggests  that  in  the  wet  case  only  a very  small  fraction 
of  the  holes  generated  in  the  Al^Oj  per  unit  dose  are  trapped  there. 

An  alternative  to  the  above  explanation  is  that  identical  traps 
and  trap  densities  exist  in  the  sapphire  for  wet  and  dry  devices. 

Data  of  Figure  45  could  then  be  explained  in  terms  of  a higher 
mobility  in  dry-oxide  devices  than  in  their  wet  counterparts.  For 
equal  trapped  charge  densities,  it  seems  reasonable  to  assume 
that  n^nv  would  be  comparable  for  the  two  device  types.  However, 
this  would  then  require  that  \x^^^{drY)  be  a factor  on  the  order  of 
10  larger  than  PbcCwet)  in  the  nonsaturation  regime  and  about  two 
orders  of  magnitude  larger  in  saturation.  This  requirement  seems 
unreasonable  and  we  judge  this  alternative  explanation  to  be  much 
less  attractive  than  the  unequal -trap -density  model  suggested  above. 

The  data  of  Figure  45  cannot  be  explained  simply  on  the  basis 
of  differing  silicon  doping  conditions.  As  mentioned  above,  the 
presence  of  an  n-type  layer  in  dry  devices  can  account  for  the  large 
pre -irradiation  value  of  compared  to  that  in  wet  units.  If  we 

assume  equal  trap  densities,  equal  trapped  charge  densities  at  a 
given  dose,  equal  back-channel  mobilities  at  a given  n-  , and  un- 
equal silicon  doping  conditions  in  the  two  cases,  an  initial  differ- 
ence in  I^  versus  dose  behavior  is  expected,  similar  to  that  shown 
in  Figure  45  at  low  doses.  However,  at  higher  doses  where  leak- 
age currents  are  significantly  larger  than  pre -irradiation  values, 
data  for  wet-oxide  devices  would  merge  with  dry-oxide  data  and  the 
saturation  value  of  would  be  the  same  in  both  cases.  Since  be- 
havior of  that  nature  is  not  observed,  we  therefore  rule  out  such 
an  explanation. 

In  order  to  determine  the  dry-to-wet  trap-density  ratio  respon- 
sible for  the  factor -of-ninety  displacement  in  Figure  45,  one  would 
need  to  relate  (or  dose)  quantitatively  to  n^^^  or  p.  The  pres- 
ent data  do  not  yield  these  relations.  It  is  interesting  to  note  that 
the  dry-to-wet  li  ratio  in  saturation  is  also  90.  However,  this  is 
fortuitous  since,  as  ahown  in  Figure  46  (devices  2 and  4),  the  satu- 
ration ratio  varies  with  temperature  due  to  the  temperature  depend- 
ence of  On  the  other  hand,  the  displacement  ratio  is  expected 

to  be  temperature  invariant  since  it  depends  only  on  the  trap-density 
ratio  under  the  present  interpretation. 
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4.  3,  3 Radiation -Induced  Reduction  of  Back -Channel  Leakage 
Current 


The  right  portion  of  Figure  44  illustrates  the  phenomenon  of 
radiation-induced  reduction  of  back-channel  leakage  current.  Re- 
ducing V^g  to  zero  after  radiation-induced  increases  in  1^^  have 
been  observed  caused  to  decrease  dramatically  to  near  its  pre- 
ir radiation  value  if  irradiation  is  continued  in  this  bias  condition. 

To  obtain  the  recovery  data  of  Figure  44  (note  break  in  dose  axis 
iust  beyond  10^  rads),  Vj^g  was  reduced  to  zero  after  10^  rads  and 
irradiation  continued.  To  measure  leakage  current  at  a given  dose, 
V^g  Wcls  momentarily  switched  to  5V  and  V to  -7V  (typically). 


This  measurement  took  2 to  3 sec  and  then  V. 


g and  V^g  were  re- 


turned to  zero.  If  the  measurement  is  not  made  rapidly,  then  I^ 
will  increase  significantly  during  the  measurement  period.  It  is 
seen  that  a few  thousand  rads  is  sufficient  to  reduce  I^  to  a satu- 
ration level. 


The  l£  reduction  process  can  be  repeated  as  evidenced  by  data 
in  Figure  47.  This  figure  shows  findings  for  a wet-oxide  transistor 
irradiated  with  V^jg  = 5V  to  a dose  of  10^  rads,  at  which  point  V^jg 
was  reduced  to  zero  and  irradiation  continued.  Substantial  recovery 
is  evident.  This  process  was  repeated  for  four  more  cycles  as 
shown.  Upon  completion  of  this  cyclying  experiment,  the  Co^® 
source  was  lowered  (dose  rate  = 0),  with  a resulting  decrease  in 
l£  to  within  25%  of  the  pre -ir radiation  value.  The  Co“  source  was 
raised  again  (with  Vds  = 0)  and  I^  increased,  and  we  attribute  this 
increase  to  a photocurrent.  (In  Figure  44,  the  final  value  of  I^ 
for  the  wet  device  following  radiation-induced  recovery  was  influ- 
enced by  the  photocurrent  and  subtracting  it  reduces  measured  Ij^ 
to  within  50%  of  its  pre -irradiation  value.)  The  recovery  phe- 
nomenon observed  here  suggests  that  periodic  reduction  of  V^jg  to 
zero  to  cause  a reduction  in  I^  is  a potential  remedy  for  leakage 
current  problems  in  SOS  devices  employed  in  a space  radiation 
environment. 


Leakage  current  studies  were  made  on  n-channel  transistors 
contained  in  CMOS/SOS  inverters  from  three  manufacturers  to  de- 
termine whether  the  radiation-induced  1^  reduction  process  ob- 
served for  custom  transistors  (Figures  44  and  47)  was  general. 
Large  variations  in  radiation  response  were  observed  from  lot  to 
lot  for  inverters.  However,  the  recovery  phenomenon  was  evident 
for  all  devices  examined.  It  is  interesting  to  note  that  in  several 


115 


cases  the  final  value  of  following  leakage  current  production  and 
subsequent  radiation -induced  reduction  was  actually  less  than  the 
pr e -ir r adiation  value.  In  one  case,  post-recovery  l£  was  an  order 
of  magnitude  smaller  than  prior  to  bombardment. 

4.3.4  Scanning -Electron -Microscope  Studies  of  Leakage  Current 
Production  and  Reduction 

Experiments  were  performed  using  a scanning  electroii  micro- 
scope (SEM)  as  a source  of  localized,  nonpenetrating  ionizing  radi- 
ation for  the  purpose  of  gaining  understanding  of  the  mechanisms  of 
leakage  current  production  and  reduction  in  irradiated  SOS  devices. 
Energy  deposition  profiles  for  low-energy  electrons  incident  on 
low-Z  materials  can  be  calculated  based  on  the  work  of  Everhart 
and  Hoff 43  and  Figures  48  and  49  show  profiles  appropriate  for  the 
geometry  of  the  present  custom  transistors  (thin  gate  metallization). 
In  Figure  48,  profiles  for  energies  from  6 to  16  keV  are  shown  and, 
for  convenience,  effective  thicknesses  of  Al,  Si,  and  AI2O3  are 
given  based  on  density  ratios  relative  to  SiO^.  It  is  seen  that  sig- 
nificant energy  deposition  occurs  in  the  sapphire  for  energies 
^9  keV.  Figure  49  shows  deposition  profiles  at  higher  energies 
and  only  the  sapphire  portion  of  a transistor  is  shown.  For  this 
figure,  distances  shown  correspond  to  actual  distances  into  the 
sapphire.  (In  discussion  below,  the  term  "electron  range"  in  the 
AI2O3  is  used  and  this  refers  here  to  that  distance  at  which 
dE/dx  = 0 in  Figure  49.) 

Irradiations  were  performed  by  raster  scanning  the  beam  over 
a device  for  an  integral  number  of  scans.  Beam  current  was  meas- 
ured by  means  of  a Faraday  cup  which  was  mechanically  moved  into 
position  between  irradiations  of  a device.  The  beam  current  was 
reduced  to  the  range  of  1 to  6 pA  by  means  of  intensity  modulation 
at  50  kHz  using  a combination  of  ac  and  dc  deflection  of  the  beam. 
Uniformity  of  irradiation  over  a device  is  assured  by  use  of  a de- 
focused  beam.  In  spite  of  the  use  of  low  beam  current  and  a de- 
focused  beam,  a significant  difference  between  Co^*^  and  SEM  ir- 
radiation conditions  is  the  much  higher  dose  rate  in  the  latter  case. 
In  the  present  study,  the  dose -rate  ratio  exceeded  10^  in  most 
cases.  Operation  was  restricted  to  a beam  current  range  where 
dose -rate  effects  on  measured  leakage  current  were  observed  to 
be  relatively  small. 

43t.E.  Everhart  and  P.H.  Hoff,  J.  Appl.  Phys . 4Z.  5837  (1971). 
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Figure  48.  Energy  deposition  profiles  for  custom  SOS  transistors 
bombarded  by  low-energy  electrons.  Profiles  were  calculated 
using  the  work  of  Everhart  and  Hoff  and  effective  thick- 
nesses of  Al,  Si,  and  AI2O3  are  shown  based  on  density  ratios 
relative  to  SiO^. 
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Figure  49.  Energy  deposition  profiles  for  the  sapphire  substrate 
of  custom  SOS  transistors  bombarded  by  low-energy  electrons. 

In  contrast  to  Figure  48,  the  x axis  corresponds  to  actual  dis- 
tances into  the  sapphire.  Note  that  the  Si-Al202  interface  is 
coincident  with  the  y axis. 

Leakage  current  studies  using  the  SEM  were  performed  on  one 
dry-oxide  chip  and  one  wet-oxide  chip,  each  of  which  contained 
several  transistors.  The  most  complete  set  of  data  were  obtained 
for  two  wet-oxide  devices  connected  in  parallel  and  irradiated  si- 
multaneously, and  we  therefore  concentrate  on  these  results  here. 
Results  for  dry-oxide  devices  were  in  qualitative  agreement  with 
wet-oxide  findings  over  the  energy  range  for  which  comparable 
data  were  obtained  (8  to  18  keV).  It  is  significant  to  note  that  the 
process  of  radiation -induced  reduction  of  leakage  current  permits 
the  performance  of  various  radiation  effects  experiments  on  a given 
device  since  the  device  can  be  returned  to  near  its  original  condi- 
tion by  simple  irradiating  with  Vds=  0.  This  procedure  was  followed 
here . 

In  one  series  of  experiments  on  wet-oxide  transistors,  the  ef- 
fectiveness of  bombarding  electrons  of  various  energies  in  giving 
rise  to  back-channel  leakage  current  was  examined.  Figure  50  shows 
leakage  current  versus  total  energy  incident  on  the  gate  electrode  for 
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Figure  50.  Leakage  current  vs  energy  incident  on  the  gate  elec- 
trode of  a wet-oxide  SOS  transistor.  Data  for  several  bom- 
barding electron  beam  energies  are  shown. 

four  beam  energies,  ranging  from  6 to  16  keV.  These  data  are 
presented  in  terms  of  total  incident  energy  since,  as  is  evident  in 
Figure  48,  nonuniform  energy  deposition  in  four  materials  is  in- 
volved and  it  is  difficult  to  establish  a more  meaningful  base  for 
comparison.  It  is  seen  that  6-keV  electrons  do  not  cause  an  increase 
in  leakage  current  whereas  the  higher -energy  cases  do.  Upon  com- 
paring the  results  of  Figure  50  with  the  profiles  of  Figures  48  and 
49,  the  following  conclusion  becomes  evident:  Energy  must  be  de- 
posited in  the  sapphire  before  radiation -induced  increases  in  back- 
channel  leakage  current  will  be  observed.  For  the  6-keV  case,  the 
energy  deposition  profile  peaks  in  the  Si02  and  the  maximum  dose 
deposited  in  the  layer  was  '^3  x 10-'  rads  (Si).  At  l6-keV,  the  peak 
of  the  profile  was  just  within  the  sapphire  and  a relatively  high 
rate  of  leakage  current  production  is  observed. 


120 


Figure  51  presents  radiation -induced  back -channel  leakage 
current  versus  dose  for  bombardment  with  30-  and  45-keV  electrons. 
For  reasons  which  are  discussed  below,  dose  in  both  cases  was  ob- 
tained using  the  average  value  of  the  energy-deposition -profile  am- 
plitude over  the  first  2 pm  of  the  sapphire.  Shown  for  comparison 
are  Co^^  data  for  a wet-oxide  transistor  taken  from  Figure  45.  The 
30-keV  data  are  in  agreement  with  45-keV  findings  and  a general 
similarity  exists  between  Co^®  and  electron  data.  Saturation  values 
of  observed  at  high  doses  are  quite  similar  and  the  rates  of  leak- 
age current  production  are  at  least  comparable.  The  results  illus- 
trated in  Figure  51  lead  us  to  a second  conclusion:  Energy  deposition 
in  the  bulk  of  the  sapphire  (i.e.,  at  distances  greater  than  several 
pm  from  the  Si-Al^O^  interface)  is  unimportant  in  terms  of  giving 
rise  to  back-channel  leakage  current.  This  point  is  clearly  illustra- 
ted by  considering  the  30-keV  profile  of  Figure  49.  Electrons  of  this 
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Figure  51.  Radiation -induced  back -channel  leakage  current  vs 

dose  for  wet-oxide  SOS  transistors.  The  Co^®  data  shown  are 
from  Figure  45.  For  electron  bombardment,  the  dose  shown 
corresponds  to  the  average  energy  deposited  in  the  first  2 pm 
of  the  sapphire  substrate. 


121 


i 


energy  only  penetrate  a few  pm  into  the  sapphire  but  give  rise  to 
a saturation  nearly  equal  to  that  produced  by  gamma  rays 

which  penetrate  through  the  entire  sapphire  thickness  (several 
hundred  pm). 

Further  insight  regarding  leakage  current  production  is  gained 
by  considering  the  electron  energy  dependence  of  saturation  leakage 
current  as  shown  in  Figure  52.  This  figure  presents  saturation 
values  as  a function  of  electron  range  ir.  the  sapphire.  (Shown  for 
comparison  is  the  mean  value  obtained  from  Co®^  measurements 
on  a set  of  transistors  from  the  same  lot.)  These  results  suggest 
a model  in  which  the  relevant  hole  traps  are  distributed  to  a depth 
of  several  pm  within  the  sapphire,  where  only  a fraction  of  these 
traps  can  be  populated  at  lower  beam  energies.  However,  the  al- 
ternative hypothesis,  in  which  the  traps  are  located  at  the  Si-Al203 
interface  (just  inside  the  AI2O3),  cannot  be  ruled  out  by  these  data 
because,  as  discussed  below,  neutralization  of  trapped  charge  oc- 
curs simultaneously  with  the  production  process.  These  two  proc- 
esses have  different  beam  energy  dependences  and  come  into  equi- 
librium in  saturation  at  each  beam  energy.  A fit  based  on  these 
competing  processes  is  shown  in  Figure  52,  which  demonstrates, 
as  will  be  seen,  that  the  observed  saturation  behavior  can  be  ex- 
plained without  invoking  spatially  distributed  traps  within  the 
sapphire , 

For  the  purpose  of  further  examining  the  region  of  importance 
in  the  sapphire,  the  effectiveness  of  various  beam  energies  in  giving 
rise  to  leakage  current  was  considered.  We  define  a relative  pro- 
duction effectiveness  (RPE)  as  AI^/E^j/Ijjdt,  where  E^^  and  are 
beam  energy  and  beam  current,  respectively,  and  is  the  change 
in  leakage  current  produced  by  a given  total  incident  energy.  Fig- 
ure 53  presents  R.PE  versus  electron  beam  energy  (open  circles) 
for  a wet-oxide  transistor.  A peak  in  effectiveness  is  observed  in 
the  range  I6  to  18  keV,  corresponding  to  beam  energies  at  which 
the  deposition  profile  peaks  at  some  slight  distance  into  the  sap- 
phire. The  sharp  decrease  in  RPE  evident  in  Figure  53  at  low  beam 
energies  is  primarily  due  to  the  smaller  fraction  of  total  incident 
energy  which  reaches  the  sapphire.  The  decrease  in  RPE  evident 
at  higher  beam  energies  is  consistent  with  either  of  the  following 
interpretations:  (1)  a finite  spatial  distribution  of  traps  which  are 
effective  in  increasing  extends  into  the  sapphire  a distance  on 
the  order  of  2 pm,  or  (2)  an  effective  generation  region  (or  collec- 
tion volume)  extends  ''-2  pm  into  the  AI2O3.  In  the  latter  case. 
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Figure  52.  Saturation  leakage  current  vs  electron  range  in  sap- 
phire for  a wet-oxide  SOS  transistor  irradiated  with  electrons 
of  various  energies.  The  mean  Co^®  saturation  value  (averaged 
over  several  samples)  is  also  shown  for  comparison.  The 
electron  ranges  shown  correspond  to  distances  at  which  dE/dx  = 0 
in  Figure  48.  The  analytical  fit  is  based  on  simultaneous  pro- 
duction and  reduction  processes  (see  text). 
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Figure  53.  Leakage -current  production  effectiveness  vs  electron 
beam  energy  for  an  irradiated  wet-oxide  SOS  transistor.  The 
fit  shown  is  based  on  simultaneous  production  and  reduction 
processes  (see  text). 
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generated  holes  would  drift  to  traps  located  at  the  Si-Al203  inter- 
face. The  important  point  to  note  is  that  the  first  few  pm  of  the 
sapphire  are  shown  here  to  be  important  in  the  leakage  current  pro- 
duction process.  It  was  on  this  basis  that  we  chose  to  express  dose 
for  the  30-  and  45-keV  data  in  Figure  51  in  terms  of  energy  deposited 
in  the  first  two  pm. 

Since  energy  deposited  in  the  sapphire  is  the  important  quantity 
in  terms  of  leakage  current  production,  it  should  be  meaningful  to 
express  RPE  in  absolute  terms  by  multiplying  the  denominator  of  the 
R PE  expression  by  that  fraction  f(E)  of  incident  energy  which  is  de- 
posited in  the  sapphire.  (The  upper  curve  in  Figure  54  shows  f(E) 
versus  electron  beam  energy.)  Absolute  production  effectiveness 
(APE),  determined  in  this  manner,  is  shown  in  Figure  53.  The  drop- 
off in  effectiveness  at  low  beam  energies  was  unexpected  since  the 
APE  should  be  independent  of  energy  at  low  beam  energies  if  f(E)  is 
the  governing  quantity.  Several  factors  could  be  responsible  for  ob- 
served behavior.  Experimental  uncertainties  are  considerably 
larger  for  the  low-energy  data  points  due  to  the  extreme  sensitivity 
of  f(E)  at  such  energies  to  assumed  silicon  thickness.  Additionally, 
production  effectiveness  had  to  be  evaluated  differently  at  d keV  than 
in  higher  energy  cases  since  a comparable  range  in  leakage  current 
was  not  covered  (i.e.,  AI^  was  different).  In  view  of  the  nonlinear 
dependence  of  on  dose,  this  implies  an  additional  uncertainty  in 
the  9-keV  case.  The  above-mentioned  competition  between  produc- 
tion and  reduction  processes  may  also  be  present,  however,  and 
would  yield  a dropoff  in  production  effectiveness  at  low  beam  energies. 
A fit  to  the  RPE  curve  based  on  such  a model  is  shown  (discussed 
further  below).  The  fit  does  not  take  into  account  the  finite  range 
of  the  production  process,  and  thus  correspondence  with  data  at 
high  beam  energies  is  not  expected. 

Studies  of  leakage  current  reduction  were  performed  as  a func- 
tion of  bombarding  electron  beam  energy,  and  Figure  54  shows  re- 
duction effectiveness  versus  energy  (open  circles)  for  a wet-oxide 
transistor.  Reduction  effectiveness  is  defined  similarly  to  produc- 
tion effectiveness  in  that  it  is  a measure  of  the  total  incident  energy 
required  to  reduce  by  a certain  (constant)  amount.  The  experi- 
mental procedure  at  a given  beam  energy  was  to  increase  by 
bombardment  with  V^jg  = 5V  and  then  to  decrease  by  zero -bias 
irradiation.  As  in  the  case  of  production  effectiveness  (Figure  53), 
reduction  effectiveness  peaks  at  an  intermediate  energy,  in  this 
case  about  11  to  12  keV.  However,  the  decrease  in  reduction  effec- 
tiveness at  higher  beam  energies  is  much  stronger  than  in  the 
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Figure  54.  Leakage -cur rent  reduction  effectiveness  (left  scale) 
vs  bombarding  electron  beam  energy  for  a wet-oxide  SOS  tran 
sistor.  A fit  to  the  data  based  on  the  amplitude  of  the  energy 
deposition  profile  at  the  Si-Al203  interface  is  shown.  Also 
shown  (upper  curve,  right  side)  is  the  fraction  of  energy  de- 
posited in  AI2O0  based  on  the  profiles  of  Figures  48  and  49. 
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production  case  which  suggests  that  energy  deposition  in  the  sap- 
phire is  not  the  governing  parameter  in  the  reduction  process.  On 
the  other  hand,  an  excellent  fit  to  the  data  is  obtained  using  the  ampli- 
tude of  the  energy  deposition  profile  at  the  Si-Al202  interface.  This 
suggests  that  electron  injection  into  the  sapphire  is  the  dominant  leak- 
age current  reduction  mechanism.  If  traps  are  distributed  within 

• the  sapphire,  that  conclusion  is  inescapable.  If,  however,  traps 
are  located  at  the  interface,  the  above  finding  does  not  rule  out 

• radiation-induced  depopulation  as  a contributing  mechanism. 

• 

Having  described  the  energy  dependence  of  the  reduction  proc- 
ess, we  now  return  to  consideration  of  the  energy  dependence  of 
saturation  leakage  current  and  relative  production  effectiveness 
(Figures  52  and  53).  Whereas  with  zero  bias  only  the  reduction 
process  is  operative,  under  bias  both  the  production  and  reduction 
processes  occur  simultaneously,  the  latter  coming  into  play  as 
saturation  is  approached.  If  the  former  is  governed  by  f(E)  and 
the  latter  by  dE/dx  (evaluated  at  the  interface),  then  both  satura- 
tion leakage  current  and  production  effectiveness  should  be  described 
by  an  appropriate  weighting  of  these  two  terms.  The  fits  in  Fig- 
ures 52  and  53  were  obtained  in  this  way,  and  they  show  the  expected 
decrease  in  saturation  leakage  current  and  production  effectiveness 
at  low  beam  energies.  Since  neither  of  these  sets  of  data  can  be  fit 
based  on  f(E)  alone  and  since  the  fits  shown  are  reasonable,  the 
model  in  which  traps  are  located  at  the  interface  appears  just  as 
tenable  as  the  one  seemingly  implied  by  the  saturation  data,  namely 
that  the  traps  are  distributed. 

An  additional  experiment  was  performed  to  examine  the  reduc- 
tion process.  Leakage  current  induced  by  means  of  a relatively 
penetrating  30-keV  beam  was  reduced  by  means  of  a 12-keV  zero- 
bias  irradiation.  Full  recovery  was  observed,  which  confirms  that 
if  hole  traps  are  distributed,  electron  injection  is  the  dominant 
recovery  mechanism.  Trap  depopulation  would  not  result  in  com- 
plete elimination  of  distributed  trapped  charge  using  the  less- 
penetrating  beam. 

4. 4 DISCUSSION 

Results  of  the  present  study  suggest  the  following  model  for 
radiation-induced  production  and  reduction  of  back-channel  leakage 
current  in  n-channel  SOS  transistors.  Ionizing  radiation  which 
penetrates  into  the  sapphire  substrate  generates  electron-hole  pairs 
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therein.  When  a drain-to -source  bias  is  applied,  a fringing  field 
exists  in  the  sapphire.  By  analogy  with  phenomena  in  irradiated 
SiO^  layers,  we  postulate  that  this  fringing  field  serves  to  sweep 
electrons  out  of  the  AI2O3  while  holes  are  trapped.  For  the  case 
of  spatially  distributed  traps,  the  holes  will  be  trapped  very  near 
their  point  of  generation.  If  all  the  traps  are  at  the  Si-Al203  inter- 
face, then  the  fringing  field  will  drive  a portion  of  the  holes  to  these 
sites  where  they  are  trapped.  The  fraction  of  generated  pairs  which 
recombines  before  electron  sweepout  occurs  should  decrease  with 
increasing  fringing  field  (i.e.,  increasing  V^jg),  and  thus  the  density 
of  trapped  holes  should  also  increase.  That  is,  should  increase 
''^th  Vjjg  in  agreement  with  experiment. As  the  dose  deposited 
in  the  sapphire  increases  for  a fixed  value  of  Vds.  will  increase 
and  reach  saturation  when  all  hole  traps  are  filled,  or,  more  gen- 
erally, when  a balance  is  reached  between  trapping  and  detrapping 
events.  If  irradiation  is  continued  to  high  doses  (^10^  rads(Si)  ), 
negatively  charged  states  will  form  at  the  Si-Al^O^  interface  which 
will  serve  to  reduce  I^  below  its  saturation  value. 

For  a given  value  of  Ai^g,  the  fringing  field  strength  decreases 
with  distance  into  the  sapphire.  Hence,  under  the  above  description, 
the  resulting  trapped  hole  concentration  would  also  decrease  with 
distance  into  the  sapphire  under  the  assumption  of  a spatial  trap 
distribution.  The  present  results  reveal  that  the  first  Z pm  of  the 
sapphire  is  the  important  region  in  the  I^  production  process.  How- 
ever, it  is  not  yet  clear  whether  this  depth  is  determined  by  a fall- 
off  in  pair  generation  probability,  as  governed  by  the  fringing  field, 
or  by  the  range  of  traps  in  the  sapphire.  (This  could  be  determined 
by  repeating  the  present  experiments  with  different  drain  biases.) 
After  leakage -current  saturation  is  achieved  by  irradiation  to  some 
high  dose  with  V^jg  applied,  continuing  the  irradiation  with  V^jg  = 0 
causes  Ij^  to  reduce  to  near  its  pre -irradiation  value.  If  the  traps 
are  spatially  distributed,  then  the  recovery  mechanism  is  injection 
of  electrons  from  the  silicon  into  the  sapphire  where  they  recombine 
with  trapped  holes.  Injected  electrons  drift  in  the  Coulomb  field 
of  the  trapped  charge  and  thus  can  neutralize  holes  trapped  some 
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distance  from  the  interface.  If  the  traps  are  located  at  the  inter- 
face, then  the  recovery  mechanisms  are  electron  injection  and/or 
radiation-induced  trap  depopulation. 

38  47  48 

Harari,  Wang  and  Royce,  and  Williams  et  al.  have  per- 
formed optical  experiments  in  which  radiation -induced  leakage  cur- 
rent in  n-channel  SOS  transistors  was  significantly  reduced  by  il- 
luminating the  sapphire  substrate  under  conditions  of  zero  drain 

bias.  Recovery  under  3.6  eV  illumination  was  ascribed  by  Harari 

47  48 

to  electron  injection  from  Si  into  Al^O^.  Wang  and  Williams 
observed  recovery  under  2.  5 eV  illumination  which  was  attributed 
to  photodepopulation  of  trapped  holes.  In  principle,  it  is  possible 
to  isolate  the  depopulation  process  from  the  injection  process  using 
optical  methods.  For  electron-beam  excitation,  both  recovery 
mechanisms  are  energetically  allowed  under  all  experimental 
conditions . 

49 

Habing  and  Shafer  observed  radiation -induced  reduction  of 
built-up  positive  charge  in  the  gate  oxide  of  CMOS -on -bulk -silicon 
devices  when  devices  previously  irradiated  xuider  positive  gate  bias 
were  further  irradiated  with  zero  gate  bias.  They  postulated  elec- 
tron injection  from  Si  into  Si02  as  the  recovery  mechanism,  and 
their  results  appear  to  be  analogous  to  the  present  observation  of 
back -channel  leakage  current  reduction. 

The  most  attractive  phenomenological  explanation  for  differ- 
ences between  the  radiation  response  of  wet-  and  dry-oxide  de- 
vices was  judged  here  to  be  the  existence  of  a larger  trap  density 
in  the  dry-oxide  case.  The  only  processing  difference  was  the  gate 
oxide  growth  temperature  and  time.  EerNisse  and  Derbenwick^O 
have  developed  a model  for  the  radiation  sensitivity  of  the  Si02-Si 
system  involving  interfacial  stress.  They  postulate  that  produc- 
tion of  hole  traps  in  the  oxide  is  related  to  viscous  shear  flow.  It 
is  not  presently  clear  whether  their  model  could  account  for 
significant  hole  trap  density  differences  in  two  Si-Al203  systems 
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which  experienced  somewhat  different  high-temperature  treat- 
ments (925°C  for  60  min  versus  1000°C  for  150  min).  This  appears 
to  be  a worthwhile  point  to  pursue  since  if  their  model  were  success- 
ful in  explaining  such  differences  it  then  might  become  clear  how 
to  optimize  SOS  processing  in  terms  of  minimizing  radiation -induced 
leakage  current. 

4. 5 CONCLUSIONS 

Major  observations  made  and  conclusions  reached  in  this  study 
are  the  following:  (a)  the  same  functional  dependence  of  radiation- 
induced  leakage  current  on  ionizing  dose  is  exhibited  for  devices 
with  wet  and  dry  gate  oxides;  (b)  differences  in  the  radiation  sensi- 
tivity of  wet-  and  dry-oxide  devices  were  attributed  to  a larger 
density  of  hole  traps  in  the  sapphire  for  dry-oxide  devices;  (c)  the 
temperature  dependence  of  leakage  current  was  attributed  to  that  of 
back-channel  carrier  mobility  in  all  cases  except  for  unirradiated 
wet-oxide  devices  at  room  temperature  and  above  where  thermal 
generation  dominates;  (d)  reduction  of  radiation-induced  leakage 
current  to  near  its  pre -irradiation  value  can  be  readily  accomplished 
by  continuing  to  irradiate  a device  with  Vds=0  ; (e)  energy  must  be 
deposited  in  the  sapphire  before  radiation-induced  increases  in  leak- 
age current  will  be  observed;  (f)  energy  deposition  in  the  sapphire 
bulk  (i.e.,  several  pm  from  the  interface)  is  unimportant  in  terms 
of  giving  rise  to  leakage  current;  the  region  which  is  effective  in 
producing  such  currents  extends  into  the  sapphire  a distance  on  the 
order  of  2 pm  from  the  interface  for  a drain  bias  of  5V:  (g)  if  hole 
traps  are  spatially  distributed,  injection  of  electrons  from  Si  into 
■^^2^3  dominant  recovery  mechanism  in  the  process  of 

radiation-induced  reduction  of  leakage  current. 
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MOS  HARDENING  APPROACHES  FOR  LOW-TEMPERATURE 

APPLICATIONS* 


5.1  INTRODUCTION 

Various  low -temperature  applications  of  MOS  devices  have 
arisen  in  the  past  few  years,  includirg  amplification,  signal  proc- 
essing, and  infrared  imaging.  Charge  transfer  devices  (CTDs) 
and  MOS  transistors  are  used  in  these  applications  at  operating 
temperatures  less  than  lOOK.  If  such  devices  are  also  employed 
in  an  ionizing  radiation  environment,  behavior  can  differ  appreci- 
ably from  that  for  room-temperature  operation.  In  particular, 
theshold  voltage  shifts  in  irradiated  MOS  transistors  (or  flatbana 
voltage  shifts  in  MOS  capacitors)  are  significantly  larger  at  low 
temperatures  than  at  room  temperature  for  devices  containing  SiO^ 
gate  insulators.  51-56  reason  for  this  is  that  holes  generated 

in  the  oxide  by  ionizing  radiation  are  relatively  immobile  at  low 
temperatures  except  when  the  electric  field  applied  to  the  oxide 
during  irradiation  is  high  (>2  x 10^  V/cm).55  At  lower  fields, 
those  electrons  generated  during  irradiation  that  escape  recombi- 
nation are  swept  out  of  the  oxide,  leaving  behind  a trapped  hole  pop- 
ulation which  gives  rise  to  large  threshold  shifts.  A hole  that  does 
not  recombine  at  low  temperatures  is  trapped  very  near  its  point  of 
generation.  It  is  also  noteworthy  that  MOS  devices  which  are 
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radiation  tolerant  at  room  temperature  exhibit  significant  radiation- 
induced  changes  at  low  temperatures.^^ 

This  section  discusses  approaches  for  reducing  the  sensitivity 
of  MOS  devices  to  low-temperature  charge  buildup.  Included  are 
analytical  and  experimental  results  relating  to  potential  hardening 
techniques.  Results  presented  herein  should  be  useful  in  the  de- 
velopment of  radiation -tolerant  MOS  devices  for  low -temperature 
applications . 

5.2  HARDENING  APPROACHES 

The  basic  effects  of  total  ionizing  dose  on  an  MOS  device  are 
charge  buildup  in  the  oxide  layer  and  interface  state  production  at 
the  SiO^-Si  interface.  Both  of  these  effects  give  rise  to  flatband 
and  threshold  voltage  shifts;  additionally,  interface  state  creation 
causes  channel  mobility  reduction  in  transistors  and  increases  in 
dark  current  and  transfer  inefficiency  in  CTDs.  As  mentioned 
above,  charge  buildup  in  SiO^  films  is  significantly  more  severe 
at  low  temperatures  than  at  room  temperature.  An  understanding 
of  the  basic  mechanisms  associated  with  this  phenomenon  leads 
us  to  suggest  the  following  approaches  to  alleviating  charge  buildup 
problems  at  low  temperatures:  (1)  carefully  select  the  field  applied 
to  the  oxide;  (2)  ion  implant  the  oxide;  (3)  use  a thin  oxide;  (4)  use 
an  insulating  film  other  than  SiO^.  This  last  item  is  not  pursued 
here.  Such  an  approach  may  prove  to  be  impractical,  but  it  should 
be  noted  that  Boesch  has  observed  AI2O3  films  on  silicon  to  be  sig- 
nificantly more  radiation  tolerant  at  low  temperatures  than  their 
Si02  counterparts.^^ 

The  first  potential  hardening  approach  --  selection  of  applied 
field  --  arose  from  the  observation  that  holes  are  mobile  at  high 
fields  in  SiO^  films  irradiated  at  low  temperatures.^^  Transporting 
holes  out  of  the  oxide  reduces  the  threshold  shift  in  an  irradiated 
transistor.  We  have  explored  this  concept  through  experiment  and 
analysis  and  findings  are  presented  below.  Regarding  the  approach 
of  hardening  by  ion  implanting  the  oxide,  previous  data^^  showed 
that  improved  radiation  tolerance  can  be  achieved  by  this  method. 

In  the  present  work,  further  analysis  of  ion  implantation  effects  is 
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given,  including  a treatment  of  the  effect  of  implantation  depth  on 
oxide  hardness.  The  third  hardening  approach  --  use  of  a thin 
oxide  --  is  beneficial  for  devices  operated  at  room  temperature.^® 
This  same  procedure  is  also  applicable  for  low -temperature  oper- 
ation because  the  total  amount  of  positive  charge  trapped  through- 
out the  oxide  decreases  with  decreasing  thickness  and  therefore 
threshold  voltage  shift  will  also  decrease.  We  treat  these  three 
topics  in  turn. 

5.3  EFFECTS  OF  APPLIED  FIELD 
5.3.1  Experimental  Results 

Experiments  were  performed  to  examine  the  applied  field  de- 
pendence of  low-temperature  charge  buildup  in  SiO^  films.  Test 
vehicles  for  this  study  were  MOS  capacitors  fabricated  by  Hughes 
Aircraft  using  procedures  developed  by  Aubuchon.^*^  These  Al-gate 
devices  contained  a dry  thermal  oxide  (910  A)  on  an  n-type  silicon 
substrate  (2  to  5 fl-cm).  Irradiations  were  performed  at  both 
room  temperature  and  77  K using  the  Co®®  source  at  Northrop  Re- 
search and  Technology  Center.  The  dose  rate  for  these  experiments 
was  40  rads(Si)/s  as  determined  by  CaF2  thermoluminescent  dosim- 
eters. Flatband  voltage  was  determined  at  several  total  doses  for 
each  device  from  in  situ  capacitance -vs -voltage  measurements 
made  with  a PAR  Model  410  CV  plotter  (1  MHz). 

5 

Figure  55  shows  flatband  voltage  shift  at  10  rads  (Si)  versus 
voltage  applied  during  irradiation  for  a group  of  identical  MOS 
capacitors.  Data  are  shown  for  77  K and  room  temperature  and 
each  data  point  was  obtained  from  measurements  on  a single  device. 
That  is,  forty  devices  were  used  to  acquire  the  data.  Figure  56 
shows  similar  data  after  a total  dose  of  10°  rads  (Si)  for  the  same 
devices  represented  in  Figure  55.  For  +10V  applied,  the  77  K 
voltage  shift  at  this  dose  is  a factor  of  7.  5 larger  than  that  at  room 
temperature.  For  an  applied  voltage  somewhat  less  than  or  con- 
siderably greater  than  +10V,  the  flatband  shift  at  77  K is  signifi- 
cantly less  than  the  peak  value. 
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Figure  55.  Flatband  voltage  shift  at  10  rads  (Si)  vs  applied  voltage 
for  MOS  capacitors  irradiated  at  77  K and  at  room  temperature 
at  a dose  rate  of  40  rads(Si)/s. 


APPLIED  VOLTAGE  IvoltsI 


Figure  56.  Flatband  voltage  shift  at  10  rads  (Si)  vs  applied  voltage 
for  MOS  capacitors  irradiated  at  77  K and  at  room  temperature 
at  a dose  rate  of  40  rads(Si)/s. 
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Negative  bias  results  for  the  present  devices  (n-type  substrate) 
should  be  indicative  of  behavior  expected  for  an  enhancement-mode 
p-channel  transistor.  For  this  application,  it  is  seen  that  the  oxide 
film  studied  is  inueed  radiation  resistant  at  room  temperature.  At 
an  applied  voltage  of  -lOV,  after  10^  rads,  the  flatband  voltage 
shift  is  quite  small  (~0.3V).  On  the  other  hand,  under  the  same 
conditions  the  flatband  shift  at  77  K is  9V,  a factor  of  thirty  larger 
than  at  room  temperature.  The  present  data  again  demonstrate 
that  SiO^  films  which  are  resistant  to  ionizing  radiation  effects  at 
room  temperature  will  experience  significant  charge  buildup  at  low 
temperatures  . 

The  rather  interesting  double -hump  structure  of  low-temperature 
data  in  Figures  55  and  56  can  be  understood  qualitatively  as  follows. 
At  low  fields  (<2  x 10^  V/  cm)  a portion  of  the  electron -hole  pairs 
generated  by  radiation  in  the  oxide  film  recombine.  Electrons  es- 
caping recombination  are  swept  out  while  re’'naining  holes  are 
"frozen  in."  As  the  field  is  increased,  a larger  ''raction  of  elec- 
trons escapes  recombination  leaving  an  incrcasiiigly  larger  amount 
of  uniformly-distributed  positive  charge.  Thus,  the  flatband  shift 
AVfb  at  a given  dose  increases  with  increasing  applied  voltage.  At 
higher  fields,  the  holes  become  mobile  which  resul  s in  less  posi- 
tive charge  in  the  oxide  and  thus  gives  rise  to  the  observed  decrease 
in  AVpjij  with  increasing  bias. 

It  is  important  to  note  that  the  high-field  data  of  Figures  55  and 
56  are  only  applicable  to  the  case  of  a relatively  low’  steady-state 
ionizing  dose  rate  and  do  not  apply  at  early  times  following  pulsed 
bombardment.  This  point  is  made  clear  in  Figure  57  where  the 
fraction  of  generated  holes  collected  at  time  t after  an  ionizing 
burst  (10-ps  pulse  of  5-k;eV  electrons)  is  plotted  versus  tempera- 
ture for  several  values  of  t.  (Information  in  Figure  57  is  based 
on  data  in  Ref.  55.)  The  data  shown  are  for  a relatively  high  value 
of  applied  field  (4  x 10^  V/cm).  At  90  K,  very  few  holes  have  been 
transported  to  the  electrodes  at  early  times  following  the  burst  and 
thus  will  be  large  for  high  doses.  At  t=  10  sec,  a significant 

amotmt  of  hole  transport  has  taken  place  so  AV^^  will  be  significantly 
less  than  at  early  times.  The  data  of  Figure  55  were  obtained  after 
a 42-min  irradiation  and  therefore  considerable  hole  transport  oc- 
curred at  high  fields  during  bombardment. 
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Figure  57.  Fraction  of  generated  holes  collected  in  an  MOS 

capacitor  at  time  t following  a pulse  of  ionizing  radiation  at 
various  temperatures.  (Information  shown  here  is  based  on 
data  in  Ref.  55. ) 

5.3.2  Analysis  of  Applied  Field  Effects 

A quantitative  framework  for  describing  the  data  in  Figures  55 
and  56  is  now  presented.  As  a starting  point,  we  assume  for  an  MOS 
capacitor  (unimplanted)  in  an  ionizing  radiation  environment  that  all 
electrons  generated  in  the  oxide  which  escape  initial  recombination  ' 
are  swept  out,  leaving  behind  trapped  holes.  Consider  an  oxide 
layer  of  thickness  t^^  in  such  a device.  If  a uniform  trapped  hole 
concentration  p exists  throughout  the  oxide,  the  resulting  flatband 
shift  is  given  by 

Initial  recombination  is  usually  taken  to  be  synonymous  with 
geminate  recombination.  In  a given  situation,  columnar  (or  "track") 
recombination  may  also  be  important  (see  Refs.  60  and  45).  We 
use  the  term  "initial"  loosely  here  to  denote  the  dominant  recom- 
bination mechanism  at  early  times  following  pair  generation. 
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If  initial  recombination  of  holes  is  neglected,  the  total  positive 

charge  density  generated  by  ionizing  radiation  can  be  calculated 

by  assuming  18  eV/pair  (translates  to  7.9  x 10^^  pairs/cm^rad 

for  Si02).  Then,  in  terms  of  ionizing  dose  y in  rads(Si),  the 

quantity  p in  coul/cm^  can  be  expressed  as  1.26 

band  voltage  shift  in  Eq.  (1)  is  then  given  by  1 . 82  x 

where  t is  in  cm. 
ox 


ads  {Si), 

X 10-°  Y 
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Hole  motion  is  now  taken  into  account.  The  following  simple 
model  was  considered  first.  At  t = 0,  electron-hole  pairs  are 
generated  uniformly  in  the  oxide.  At  t = 0'*',  a portion  of  those 
pairs  have  already  recombined  and  all  of  the  remaining  electrons 
are  swept  out,  leaving  a uniform  "frozen-in"  hole  distribution. 

/ C p 

For  high  applied  fields  (>2  x 10°  V/cm),  hole  transport  occurs. 

At  such  fields,  initial  recombination  is  negligible  (discussed  be- 
low) and  thus  all  generated  holes  remain  in  the  oxide  at  t = 0'*’.  In 
this  simple  model,  the  entire  distribution  of  holes  is  assumed  to 
move  uniformly  through  the  oxide.  For  example,  under  positive 
bias  the  hole  population  moves  a "transport  distance"  d^  toward 
the  SiO^-Si  interface  in  a time  t^. 

The  model  described  above  was  pursued  and  found  to  be  un- 
satisfactory in  accounting  for  the  data  of  Figures  55  and  56.  Specif- 
ically, in  determining  those  values  of  transport  distance  d^  re- 
quired to  fit  AV^  data,  it  was  found  that  both  the  magnitude  of  d^ 
at  a given  field  and  its  functional  dependence  on  applied  field  dis- 
agreed significantly  when  results  for  positive  and  negative  bias 
were  compared.  This  physically  unrealistic  result  led  us  to  em- 
ploy a more  meaningful  model  in  which  a distribution  of  hole  mo- 
bilities is  taken  into  account.  That  is,  at  a given  time  following 
carrier  excitation  the  more  mobiles  holes  have  been  swept  out  of 
the  oxide  while  their  less  mobile  counterparts  remain  there  until 
some  later  time.  Such  a viewpoint  is  consistent  with  both  the 
multiple -trapping  model^5,6l  CTRW  model^^'  of  charge 

transport  in  SiO;?. 
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Figure  58  illustrates  this  more  realistic  transport  model. 
Shown  is  an  MOS  capacitor  with  a charge  density  p,  dependent  on 
distance  x,  superimposed  on  the  oxide.  The  solid  line  depicts  the 
charge  distribution  at  a particular  time  after  excitation  and  we 
characterize  this  distribution  in  terms  of  a transport  distance  d^. 
as  in  the  simple  model.  However,  d^  is  now  defined  differently  as 
shown  in  the  figure.  Additionally,  it  is  assumed  that  when  d^  = 
then  p is  equal  to  its  maximum  value  Pq  only  at  the  SiO^-Si  inter- 
face. That  is,  the  "slope"  of  the  charge  distribution  is  assumed 
to  remain  constant  as  charge  is  swept  out  of  the  oxide.  The  solid 
line  in  Figure  58  corresponds  to  a value  for  d less  than 

Also  shown  for  illustration  (dashed)  are  two  other  cases;  d,  = t /Z 
, , t ox 

> ‘ox/^- 

Expressions  for  flatband  voltage  shift  in  terms  of  transport 
distance  for  the  model  depicted  in  Figure  58  are  now  given.  Under 
an  applied  positive  bias,  for  d^  ^ ^ox^^' 
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Figure  58.  Illustration  of  a more  realistic  model  for  hole  motion 
at  high  fields  and  low  temperatures  in  SiO^  (positive  bias  case). 
Shown  are  three  hole  distributions  corresponding  to  three  values 
of  transport  distance  d^. 
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For  the  case  of  an  applied  negative  bias,  the  hole  distributions  are 

of  the  same  form  as  those  in  Figure  58  except  that  now  we  define 

d,  in  terms  of  distance  from  the  SiO^-Si  interface.  For  d^  s t /2, 
t , . 2 t ox 

we  obtain 


Equations  (2)  through  (5)  can  be  viewed  as  yielding  the  temporal 
dependence  of  AY^  since  d^  is  a function  of  time  following  bom- 
bardment. For  a fixed  time  following  bombardment,  variations  in 
d^  can  be  obtained  by  applying  different  values  of  field  to  a set  of 
samples  irradiated  eind  measured  under  otherwise  identical  conditions. 

55 

As  mentioned  previously,  a complication  arises  when  com- 
paring expected  and  measured  flatband  voltage  shifts  at  low  tem- 
peratures for  doses  greater  than  ~10^  rads  (Si).  At  such  doses,  the 
oxide  field  will  be  perturbed  due  to  the  substantial  charge  buildup 
that  exists.  This  perturbation  is  manifested  as  a decrease  in  field 
near  the  metal -Si02  interface  and  an  increase  in  field  near  the  Si02“ 

Si  interface.  In  the  increased  field  region,  hole  transport  will  occur 
if  the  field  there  is  greater  than  2x10  V/cm.  In  the  decreased 
field  region,  more  initial  recombination  of  generated  electron -hole 
pairs  will  occur.  Flatband  voltage  shifts  will  be  reduced  by  both  of 
these  mechaidsms  relative  to  values  expected  on  the  basis  of  a uni- 
form field  throughout  the  oxide.  Hence,  in  analysis  of  the  present 
55 

J.R.  Srour,  S.  Othmer,  O.L..  Curtis,  Jr.,  andK.Y.  Chiu, 
IEEE  Trans.  Nucl.  Sci.  23,  1513  (1976). 
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experimental  results  we  only  treat  data  obtained  at  10  rads  (Si) 
(Figure  55)  in  an  attempt  to  avoid  these  complications.  (Calcula- 
tions by  Boesch  and  McGarrity^^  for  an  oxide  of  thickness  simi- 
lar to  that  utilized  here  indicate  non-negligible  field  modification 
at  10^  rads  (Si)  for  a device  irradiated  at  80  K with  a lOV  bias. 
Assuming  that  their  calculations  apply  quantitatively  to  the  present 
experimental  conditions,  it  still  appears  likely  that  the  predicted 
field  perturbation  for  this  dose  would  not  significantly  alter  con- 
clusions made  below  which  are  based  on  data  fitting  with  the  model 
of  Figure  58. ) 

Analysis  of  low-temperature  flatband  shift  versus  applied  volt- 
age data  in  Figure  55  was  performed  using  Eqs . (1)  through  (5). 

At  10^  rads(Si)  for  an  oxide  thickness  of  910  A,  Eq.  (1)  yields  = 

15.2V.  That  is,  the  maximum  predicted  value  of  is  15.2V 

and  this  value  should  result  when  all  holes  generated  in  the  oxide 
at  10^  rads  are  trapped  without  being  transported  and  all  generated 
electrons  are  swept  out.  Measured  AV^  arises  from  holes  trap- 
ped in  the  oxide  and  the  amoimt  and  distribution  of  these  trapped 
charges  can  be  attributed  to  two  field -dependent  effects.  First,  the 
number  of  holes  available  for  contributing  to  AV^  depends  on  the 
amount  of  initial  recombination  that  has  occurred.  That  is,  the 
dependence  of  electron  yield  on  applied  field  is  important.  As  yield 
increases,  AVf^  increases.  Second,  at  high  fields  holes  are  trans- 
ported out  of  Idle  oxide,  thereby  reducing  AV£j^. 

56 

Boesch  and  McGarrity  determined  the  field  dependence  of 
electron  yield'"  for  MOS  capacitors  irradiated  at  80  K by  measuring 
AVfb  at  early  times  following  pulsed  excitation  (4-p.s  pulse  of  13- 
MeV  electrons).  For  fields  > 2 x 10^  V/cm,  they  observed  a nearly 
saturated  yield  corresponding  to  sweepout  of  nearly  all  the  elec- 
trons generated.  At  such  fields,  apparently  all  generated  electrons 
escape  initial  recombination  and  are  swept  out  of  the  oxide.  For 

“"X“" 

"For  MOS  capacitors  irradiated  at  low  temperatures,  we  prefer 
to  ascribe  the  increase  in  AV^j^  with  field  to  an  increase  in  "electron 
yield"  since  the  number  of  electrons  escaping  recombination  and  be- 
ing swept  out  of  the  oxide  increases  with  field.  One  could  also  de- 
scribe the  same  situation  in  terms  of  "hole  yield"  since  the  number 
of  trapped  holes  also  increases  with  field  (except  at  high  fields)  and 
gives  rise  to  increases  in  AV^j^.  Boesch used  the  latter 
nomenclature . 

E.  Boesch,  Jr.,  andJ.M.  McGarrity,  IEEE  Trans.  Nucl. 
Sci.  1520  (1976). 
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Figure  59.  Transport  distance  vs  applied  field  based  on  calcula- 
tions using  the  model  illustrated  in  Figure  58. 

this  regime,  AV^  can  be  attributed  to  the  entire  generated  hole 
population  modified  by  transport  according  to  the  model  of  Figure  58. 

An  iterative  procedure  was  used  in  conjunction  with  Eqs.  (2) 
and  (5)  to  calculate  transport  distance  d^  based  on  measured  values 
of  AV^j^  at  10  rads  (77  K data  of  Figure  55).  Figure  59  shows  results 
of  these  calculations  in  terms  of  dj.  versus  oxide  field.  For  fields 
>3x10  V/cm,  positive  and  negative  bias  resiilts  are  in  good  agree- 
ment. As  mentioned  above,  such  agreement  is  expected  for  a physi- 
cally meaningful  model.  These  high-field  results  are  fit  reasonably 
well  with  a unity-slope  line.  At  lower  fields,  deviation  of  the  cal- 
culated points  from  this  line  is  due  to  the  influence  of  electron  yield. 
That  is,  Eqs.  (2)  through  (5)  assume  that  all  generated  electrons 
have  been  swept  out  which  is  not  the  case  at  lower  fields.  The  unity- 
slope  fit  yields  dj  = kE,  where  k = 1.5x  10"^^  cm^V^.  This  ex- 
pression was  then  used  to  fit  positive  bias  flatband  shift  data.  Equa- 
tions (2)  and  (3)  were  utilized  to  obtain  an  appropriate  charge  trans- 
port curve  and  the  result  is  shown  in  Figure  60.  The  next  step  was 
to  empirically  fit  a yield -vs -field  curve  to  the  low-field  data.  The 
criterion  employed  was  to  obtain  a curve  that,  when  considered 
simultaneously  with  the  charge  transport  curve,  would  yield  a good 
composite  fit  to  the  experimental  data.  The  resulting  curve  is  shown 
in  Figure  60  along  with  the  composite  fit  (solid  curve).  (To  obtain 
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Figure  60.  Flatband  voltage  shift  at  77  K vs  oxide  field  for  the  case 
of  an  applied  positive  bias.  Results  of  data  fitting  using  the 
charge  transport  model  of  Figure  58  in  conjunction  with  a 
yield-vs -field  component  are  shown. 


this  fit,  the  charge  transport  curve  is  multiplied  by  the  yield-vs - 
field  curve  normalized  to  its  saturation  value.)  By  considering  the 
charge -transport  and  the  yield -vs -field  mechanisms  jointly,  an 
excellent  fit  to  the  applied  field  dependence  of  flatband  voltage  shift 
at  77  K is  obtained. 


The  same  fitting  procedure  was  also  applied  to  negative  bias 
data.  A charge  transport  curve  fit  the  highest-field  data  but  de- 
viated from  the  data  at  intermediate  fields.  Reducing  k to  1.4  x 
10 12  improved  the  transport  fit  somewhat.  A rather  severe  linear 
yield-vs -field  curve  was  required,  in  conjunction  with  either  of 
these  transport  curves,  to  obtain  a fit  for  data  at  fields  < 2 x 
10^  V/cm.  A composite  fit  exhibited  a good  fit  to  the  data  for  fields 
^ 10  V/cm  and  > 3 x 10  V/cm.  The  peak  of  the  composite  fit  oc- 
curred at  a somewhat  lower  field  than  exhibited  by  the  data.  Al- 
though the  general  features  of  the  negative  bias  data  are  well  de- 
scribed by  the  modeling  approach  employed,  the  fit  obtained  was  not 
as  good  as  that  for  positive  bias. 
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The  influence  of  an  initial  transport  distance  d,,  on  the  data 
fitting  procedure  was  considered.  Boesch  and  McGarrity^“  observed 
that  initial  flatband  shifts  at  80  K differed  slightly  depending  on  the 
polarity  of  the  applied  bias.  They  account  for  this  observation  in 
terms  of  an  initial  transport  of  holes  in  the  valence  band  subsequent 
to  thermalization  but  before  the  first  capture  event.  It  seems  rea- 
sonable to  expect  that  d^  will  vary  directly  with  applied  field  and 
thus  correspond  to  an  initial  Schubweg.  Treated  in  this  manner,  d^ 
was  found  to  have  negligible  effect  on  the  fitting  procedures  used  here 
for  positive  bias  data.  For  negative  bias,  an  initial  Schubweg  has 
an  important  effect.  However,  such  a modification  to  the  fitting  pro- 
cedure actually  caused  the  fit  at  high  fields  to  be  worse  than  without 
this  correction. 


It  is  of  interest  to  compare  the  yield -vs -field  curve  in  Figure  60 
with  previous  findings.  Figure  6l  shows  this  curve  along  with  that 
obtained  by  Boesch.  The  shape  of  the  two  curves  is  nearly  the 
same  and  excellent  agreement  is  noted  for  fields  ^ 2 x 10^  V/cm. 

This  agreement  supports  the  validity  of  the  present  modeling  approach. 
Boesch  obtained  the  solid  curve  of  Figure  6lby  measuring 
shortly  ('^-l  ms)  after  a pulse  of  ionizing  radiation.  At  such  times, 
the  dependence  of  AVfb  on  applied  field  follows  the  electron  yield -vs - 
field  curve  because  hole  transport  is  negligible.  At  much  later  times, 
or  in  a steady-state  ionizing  radiation  environment,  data  such  as  that 
in  Figures  55  and  56  are  obtained. 


It  is  perhaps  surprising  that  transport  distance  was  observed  to 
depend  linearly  on  field  for  data  of  both  polarities  in  Figure  59.  Such 
a dependence  yields  an  effective  hole  mobility  that  is  independent  of 
field.  To  obtain  a rough  estimate  of  mobility,  a characteristic  time 
of  21  min.  (i.e.,  one-half  the  time  it  takes  to  deliver  10^  rads  at  a 
rate  of  40  rads/s)  was  assumed.  From  Figure  59,  d(dj.)/dE  = = 

1.5  X 10"^^  cm^V^,  which  yields  an  estimated  mobility  of 
~10"^^  cm^/Vsec.  This  observation  is  in  contrast  to  previous  re- 
sults in  which  mobility  was  noted  to  depend  on  at  intermedi- 

ate temperatures  and  on  E^^  at  low  temperature  (90K).  However, 


55j.r.  Srour,  S.  Othmer,  O.L.  Curtis,  Jr.,  and  K.Y.  Chiu, 
IEEE  Trans.  Nucl.  Sci.  1513  (1976). 

56h.E.  Boesch,  Jr.,  and  J.M.  McGarrity,  IEEE  Trans.  Nucl. 
Sci.  Z3.  1520  (1976). 

61  O.L.  Curtis,  Jr.,  and  J.  R. Srour , J.  Appl.  Phys . 48,  3819 
(1977). 
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Figure  61.  Fractional  yield  vs  oxide  field  for  the  present  study 
(from  Figure  60)  and  for  that  of  Boesch  and  McGarrity.^^ 

mobility  was  derived  in  a different  manner  than  in  the  cited  work. 

In  the  present  case,  a quantity  related  to  carrier  transport  was 
measured  at  a fixed  time  following  bombardment  whereas  in  pre- 
vious determinations  a fixed  amount  of  transported  charge  was  used 
as  the  basis  for  evaluating  mobility.  We  have  found  that  evaluating 
field -dependent  transport  data  by  these  two  methods  can  yield  sig- 
nificantly different  results  and  this  feature  may  be  responsible  for 
observations  of  a field -independent  mobility  in  one  case  and  a strong 
field  dependence  in  other  cases.  The  fact  that  test  specimens  used 
in  the  present  study  were  different  from  those  employed  in  previous 
work  may  also  be  important.  The  most  important  observation  to 
be  made,  however,  is  that  the  region  of  apparent  linearity  in  Fig- 
ure 59  is  fairly  narrow.  At  low  fields  the  dependence  of  carrier 
yield  on  applied  field  is  imporant  and  at  high  fields  saturation  of  d^ 
is  expected.  These  two  effects  possibly  mask  a nonlinear  depend- 
ence of  transport  distance  on  field  and,  thus,  a field -dependent 
mobility . 

56h.E.  Boesch,  Jr.,  and  J.M.  McGarrity,  IEEE  Trans.  Nucl. 
Sci.  2^,  1520  (1976). 
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5.4  EFFECTS  OF  ION  IMPLANTATION 


A reasonable  approach  to  reducing  the  effects  of  positive 
charge  buildup  in  SiO^  films  irradiated  at  low  temperatures  is  to 
introduce  electron  traps  for  compensation.  Electron  (and  hole) 
traps  can  be  introduced  into  SiO^  by  ion  implantation^^  and  we  have 
previously  examined  charge  buildup  at  77  K in  an  Al'*’-implanted 
specimen.  A factor  of  two  improvement  was  observed  compared 
to  an  unimplanted  device.  The  implantation  employed  was  not  opti- 
mum since  traps  were  only  introduced  into  that  half  of  the  oxide 
nearer  the  Al-SiO^  interface.  Analysis  of  ion  implantation  effects 
is  presented  here  which  is  useful  for  defining  optimum  structures 
and  for  predicting  low-temperature  radiation  response.  (This 
analysis  is  an  extension  of  previous  work^^to  low-temperature 
conditions  . ) 

Consider  an  MOS  capacitor  with  an  ion-implanted  oxide.  We 
denote  the  implantation  depth  as  dj^  (measured  from  the  metal -SiO^ 
interface)  and  the  region  between  this  interface  and  d^  is  assumed 
to  be  uniformly  damaged,  i.e.,  to  contain  a uniform  trap  concen- 
tration. The  mean  electron  range  (Schubweg)  in  the  implanted 
zone  between  the  metal-Si02  interface  and  dj^  is  given  by  E, 

where  and  are  mobility  and  lifetime,  respectively,  and  E is 
applied  electric  field.  It  is  assumed  that  at  low  temperatures  all 
holes  generated  in  the  oxide  by  ionizing  radiation  in  both  the  im- 
planted and  unimplanted  zones  are  immobile.  All  electrons  gen- 
erated in  the  unimplanted  zone  are  assumed  to  be  completely  swept 
out  of  that  region,  either  into  the  implanted  zone  or  into  the  silicon 
depending  on  the  polarity  of  the  applied  bias.  The  case  of  a posi- 
tive bias  applied  to  the  metal  electrode  is  considered  here. 

Electrons  generated  in  the  implanted  region  will  move  an  aver- 
age distance  lUg  toward  the  metal  electrode  before  being  perman- 
ently trapped.  Thus,  the  average  position  of  these  trapped  elec- 
trons is  a distance  (d^  - U)g)/2  from  the  metal -SiO^  interface. 
Electrons  generated  in  the  unimplanted  region  of  the  oxide  will  be 
swept  toward  the  implanted  zone  and  move  an  average  distance 
into  this  region  before  being  permanently  trapped.  The  average 
position  of  these  trapped  electrons  is  a distance  d^  - tUg  from  the 

^^J.R.  Srour,  S.  Othmer,  O.L.  Curtis,  Jr.,  andK.Y.  Chiu, 
IEEE  Trans.  Nucl.  Sci.  1513  (1976). 

^^O.L.  Curtis,  Jr.,  J.R,  Srour,  and  K.Y.  Chiu,  IEEE  Trans. 
Nucl.  Sci.  2^,  2174  (1975). 
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metal-Si02  interface.  Thus,  electrons  trapped  at  two  average 
positions  in  the  iniplanted  region  will  compensate  for  a portion  of 
the  positive  charge  which,  by  assumption,  is  trapped  uniformly 
throughout  the  oxide  . 

It  is  of  interest  to  compare  expected  flatband  voltage  shifts 
for  implanted  and  unimplanted  MOS  capacitors.  For  an  unim- 
planted unit,  the  average  position  of  the  trapped  positive  charge 
is  at  For  an  implanted  sample,  the  two  contributions  to 

AVfjj  from  trapped  electrons  are  subtracted  from  the  trapped  hole 
contribution.  In  addition,  different  fractions  of  the  total  amount 
of  negative  charge  generated  are  trapped  at  the  two  average  posi- 
tions. It  can  be  shown  that  the  ratio  of  flatband  voltage  shifts  for 
implanted  and  unimplanted  specimens  is  given  by 


AV_  (implanted) 
fb 

AV^  (unimplanted ) 


= 1 - 


d.  - UJ 
i e 


) -( 


d.  + UJ 


Equation  (6)  relates  the  expected  reduction  in  AV^^  accom- 
plished by  implantation  to  both  the  Schubweg  in  the  implanted  zone 
and  the  implantation  depth.  In  order  to  perform  example  calcula- 
tions, a value  for  uUg  is  needed.  Previous  work^^  yielded  a pg  T 
value  at  room  temperature  of  "'-2  x 10“^^  cm^V”^  for  the  implanted 
zone  of  20-keV  A1  -implanted  specimens  (fluence  = 10^^  cm"^). 

For  an  applied  field  of  10^  V/cm,  this  value  translates  to  a value 
for  iDg  of  200  A.  We  assume  that  UJ^  is  temperature  independent 
and  thus  that  the  value  obtained  at  room  temperature  can  be 

applied  at  low  temperatures.  The  primary  temperature  depend- 
ence expected  is  that  of  capture  cross-section  and  this  considera- 
tion is  neglected  here  for  simplicity.  Figure  62  presents  results 
of  calculations  performed  using  Eq.  (6).  The  quantity  AV^(im- 
planted / AV£^ (unimplanted)  is  plotted  versus  implantation  depth  for 
three  values  of  u)  and  an  oxide  thickness  of  1000  A.  The  U)  value 

e C 0 

of  ZOO  A was  chosen  to  be  representative  of  conditions  for  a field 
cf  10®  V/cm  and  an  implantation  fluence  of  10^^  cm”^.  The  lower 
values  of  were  selected  to  examine  the  expected  effects  of  higher 
fluences  (i.e.,  shorter  electron  Schubweg).  The  intuitively  ob- 
vious results  were  obtained:  increasing  the  ratio  of  implantation 
depth  to  oxide  thickness  and/or  reducing  the  electron  Schubweg 

^^O.L.  Curtis,  Jr.,  J.R.  Srour,  andK.Y.  Chiu,  IEEE  Trans. 
Nucl.  Sci.  2174  (1975). 


reduces  for  an  implanted  specimen  relative  to  that  for  an  un- 

implanted unit.  For  example,  implanting  ions  to  a depth  of  800  A 
at  such  a fluence  to  cause  U)g  to  equal  ZO  A is  predicted  to  reduce 
AV£j^  by  a factor  of  13. 

As  the  applied  field  is  increased,  will  increase  with  the 
result  that  a larger  fraction  of  the  electrons  generated  escape  being 
trapped  in  the  implanted  zone.  This  effect  will  diminish  the  degree 
of  compensation  of  trapped  positive  charge.  However,  as  discussed 
above,  holes  are  mobile  at  high  fields  in  unimplanted  SiOp  films  at 
low  temperatures.  In  a more  complete  analysis  of  charge  compen- 
sation and  charge  buildup  in  implanted  specimens  at  low  tempera- 
tures, one  should  consider  motion  of  holes  out  of  the  unimplanted 
zone  into  either  the  silicon  substrate  or  into  the  implanted  region, 
depending  on  the  polarity  of  the  applied  bias.  In  addition  to  neglec- 
ting hole  transport,  initial  recombination  of  generated  electrons 
and  holes  was  not  considered  in  the  above  discussion.  The  assump- 
tion was  made  in  previous  work^^that  ion  implantation  has  no  ef- 
fect on  the  processes  of  generation  and  initial  recombination  of 
electron -hole  pairs,  and  that  same  assumption  was  made  here  in 
developing  Eq.  (6). 

5.5  EFFECTS  OF  OXIDE  THICKNESS 

The  degree  of  radiation  tolerance  achievable  at  low  tempera- 
tures by  reducing  the  oxide  thickness  is  now  considered.  Decreasing 
t^jj.  increases  the  oxide  field  for  a constant  applied  voltage.  Thus, 
AV^b  is  reduced  in  two  ways:  (1)  the  total  number  of  holes  generated 
in  the  oxide  decreases;  (2)  charge  transport  increases.  Sander  and 
Gregory  observed  threshold  voltage  shift  to  vary  as  t^^  for  a set 
of  MOS  transistors  irradiated  at  76  K.  The  model  they  used  to  ac- 
count for  this  observation  is  that  which  leads  to  Eq.  (1).  The  pres- 
ent work  reveals  that  for  high  applied  fields  this  simple  model  is 
inadequate  since  significant  transport  of  holes  out  of  the  oxide  oc- 
curs. Equation  (1)  gives  the  maximum  flatband  shift  (AV^)  that 
can  occur  at  a given  dose  for  a fixed  oxide  thickness  and  this  max- 
imum results  when  all  electrons  generated  are  swept  out  and  no 
recombination  occurs.  Equations  (2)  through  (5)  can  also  be  re- 
written in  terms  of  AV^  . Additionally,  transport  distance  d^  is 

53h.H.  Sander  and  B.L.  Gregory,  IEEE  Trans.  Nucl.  Sci.  22, 
2157  (1975). 

^2o.L.  Curtis,  Jr.,  J.R.  Srour,  andK.Y.  Chiu,  IEEE  Trans. 
Nucl.  Sci.  ZZ,  2174  (1975). 
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given  by  kE  (Figure  59).  Also,  E = where  V is  applied 

voltage.  Thus,  Eqs . (2)  and  (3)  can  be  rewritten 

as 


OX 


ford^t  /Z(ort  ^ (2kV)^^^),  and 
t ox  ox 


AV  = AY™  (i  - 4 Ml!  + 2 0^^ 

fb  fb  I 3 4 3.6 


ox 


ox 


(7) 


(8) 


, ,1/2 

for  dj.  S t^^/Z  (or  t^^  s (2kV)  ).  These  expressions  contain 

both  the  thickness  and  applied-field  dependences  of  low-temperature 
flatband  shift. 


When  applying  Eqs.  (7)  and  (8),  it  should  be  noted  that  all  holes 
have  been  transported  out  of  the  oxide  when  d^.  = t^^.  This  condi- 
tion can  also  be  expressed  as  t^^^  = (kV)^^^,  thereby  giving  a min- 
imum value  of  t^^  for  which  (7)  and  (8)  apply.  Figure  63  shows 
results  of  calculations  using  these  equations  and  Eq.  (1).  The  solid 
curve  corresponds  to  the  maximum  shift  attainable  for  oxide  thick- 
nesses from  300  to  900  A after  10^  rads  (Si).  The  dashed  curves 
correspond  to  various  values  of  applied  voltage  and  these  curves 
reveal  the  reduction  in  that  occurs  due  to  hole  transport. 

An  example  will  illustrate  the  effectiveness  of  reducing  tQ^ 
in  increasing  radiation  tolerance.  Consider  calculations  for  900  A 
and  450  A oxides  in  Figure  63.  The  maximum  shifts  are  14.7  and 
3.7V  for  these  cases,  respectively,  and  this  reduction  is  due  to 
the  tgx  dependence  of  AV£jj.  The  reduction  in  charge  buildup  is 
much  more  dranatic  when  a field  is  applied.  For  lOV  applied, 
the  predicted  reduction  in  AV^  due  to  hole  transport  is  insignii- 
cant  (14.7V  reduces  to  14.5V)  at  900  A.  However,  at  450  i,  the 
maximum  shift  of  3.7V  is  predicted  to  be  reduced  to  0.82V.  Thus, 
the  present  calculations  indicate  that  for  an  applied  voltage  of  lOV, 
reducing  t^^  from  900  to  450  k will  decrease  AV£]^  by  a factor  of 
''-18.  If  reliable  thin  (^500  A)  oxides  can  be  fabricated  for  use  at 
applied  voltages  ^lOV,  then  it  appears  that  low-temperature 
charge  buildup  problems  will  be  largely  eliminated  in  a steady- state 
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Figure  63.  Calculated  flatband  voltage  shift  at  10  rads  (Si)  vs 
oxide  thickness  for  four  values  of  applied  voltage.  Also 
shown  (solid  curve)  is  the  maximum  flatband  shift  which 
results  when  there  is  no  recombination  and  no  hole 
transport. 


ionizing  radiation  environment.  However,  this  hardening  approach 
will  not  be  nearly  as  effective  at  early  times  following  pulsed  bom- 
bardment. For  that  case,  the  reduction  in  AV^  that  results  when 
t^^  is  reduced  will  follow  the  solid  curve  of  Figure  63. 

5.6  SUMMARY  AND  CONCLUSIONS 

We  have  found  that  the  applied  field  dependence  of  flatband 
voltage  shift  in  MOS  capacitors  irradiated  at  77  K can  be  accounted 
for  in  terms  of  the  field  dependence  of  electron  yield  and  the  trans- 
port of  holes  at  high  fields.  Analysis  of  ion  implantation  effects 
indicates  that  a significant  improvement  in  radiation  tolerance 
should  be  achievable  at  low  temperatures  by  this  method.  A simul- 
taneous consideration  of  the  effects  of  oxide  thickness  and  applied 
field  on  charge  buildup  in  an  unimplanted  oxide  suggests  that  re- 
ducing tg^  to  500  A will  largely  eliminate  low-temperature  prob- 
lems in  a steady-state  ionizing  radiation  environment  as  long  as 
the  applied  voltage  is  ^lOV.  If  such  conditions  prove  to  be  im- 
practical, then  ion  implantation  may  be  a valuable  hardening  pro- 
cedure. Use  of  a thin  oxide  with  a relatively  high  applied  field 
appears  to  be  a more  attractive  approach  to  realizing  hardened 
MOS  devices  for  low -temperature  applications  in  a steady-state 
environment.  At  early  times  following  pulsed  bombardment,  this 
approach  will  not  be  as  effective. 


SECTION  6.0 


DETERMINATION  OF  THE  DIFFUSION  LENGTH  AND 
DRIFT  MOBILITY  IN  SILICON  BY  USE  OF  A 
MODULATED  SEM  BEAM* 


6.1  INTRODUCTION 

The  parameters  of  minority-carrier  diffusion  length  and  drift 
mobility  are  of  importance  in  all  semiconductor  devices  relying  on 
minority-carrier  conduction,  such  as  bipolar  transistors  and  solar 
cells.  When  minority  carriers  are  locally  injected  into  a semi- 
conductor, the  spatial  distribution  is  governed  by  the  diffusion 
length  and  the  temporal  decay  by  the  lifetime.  These  quantitities 
are  related  by  the  drift  mobility.  The  lifetime  in  silicon  is  found 
to  be  extremely  sensitive  to  material  imperfections  and  the  con- 
centration of  certain  impurities.  For  example,  these  factors  may 
limit  the  efficiency  of  solar  cells  in  that  the  effective  collection 
volume  for  generated  carriers  is  given  by  the  diffusion  length. 
Currently,  efforts  are  underway  to  develop  technologies  for  low- 
cost  silicon  solar  cells  for  terrestrial  applications.  The  large 
cost  reduction  required  compels  the  use  of  silicon  which  is  con- 
siderably less  ideal  than  that  employed  in  space  cells.  The  mate- 
rial may  be  marginal  in  diffusion  length,  nonhomogeneous  or  even 
polycrystalline.  In  such  material,  the  drift  mobility  may  differ 
from  that  expected  in  high-grade  silicon  and  therefore  needs  to  be 
characterized  in  addition  to  lifetime  or  diffusion  length.  If  life- 
time is  measured,  moreover,  the  drift  mobility  is  required  to 
determine  the  equivalent  diffusion  length.  Since  drift  mobilities 
are  difficult  to  obtain  experimentally,  particularly  in  materials  of 
short  diffusion  length.  Hall  measurements  of  mobility  are  gener- 
ally relied  upon.  The  measured  quantity  here  is  the  majority- 
carrier  mobility.  Moreover,  such  measurements  lead  to  erron- 
eous conclusions  in  the  presence  of  the  type  of  inhomogeneities 
expected  to  characterize  low-grade  silicon.  A direct  determination 


* 

The  original  research  upon  which  this  section  is  based  was 
performed  under  previous  HDL  contracts  (DAAG39-73-C-0171  and 
DAAG39-69-C-0039).  The  review  of  that  work  given  here  was  pre 
pared  during  the  present  contract  period  and  presented  at  the  1978 
Scanning  Electron  Microscopy  Symposium. 
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of  minority -car rier  drift  mobility  would  therefore  be  of  interest. 

In  the  technique  discussed  here,  both  diffusion  length  and  drift 
mobility  may  be  determined. 

Another  problem  to  which  the  present  methods  are  applicable 
is  the  degradation  of  bipolar  transistors  under  energetic  neutron 
bombardment.  In  the  case  of  microwave  (narrow  basewidth)  tran- 
sistors, for  example,  decrease  in  drift  mobility  with  fluence  may 
be  determinative.  This  occurs  at  such  high  fluences  that  the  dif- 
fusion length  will  have  degraded  to  the  order  of  a few  pm,  requiring 
the  determination  of  transport  properties  within  spatial  resolution 
of  this  order . 

The  conventional  method  of  measuring  the  drift  mobility  is  the 
Haynes -Shockley  technique  in  which  the  transport  time  of  a pulse 
of  minority  carriers  is  monitored  under  application  of  a known 
field.  The  degradation  of  drift  mobility  in  germanium  under  neu- 
tron bombardment  has  been  studied  in  this  manner. The  method 
is  difficult  to  implement  for  material  of  short  diffusion  length, 
however,  and  the  alternative  approach  proposed  here,  which  may 
be  regarded  as  a cw -equivalent  of  the  above  pulse  technique,  is 
more  suitable  to  this  case.  A modulated  SEM  beam  is  employed 
as  a source  of  minority  carriers  and  a planar  junction  as  a collec- 
tor of  these  carriers.  Injected  carriers  are  swept  toward  or  away 
from  the  junction  by  means  of  an  applied  field  and  changes  in  signal- 
delay  per  unit  field  are  used  to  derive  the  mobility. 

Since  the  above  technique  requires  special  sample  geometries 
and  unusual  experimental  facilities,  it  is  not  likely  to  find  wide 
application.  Comparison  of  the  above  method  was  therefore  made 
with  field-free  methods  of  determining  the  drift  mobility  which 
could  readily  be  performed  elsewhere.  These  methods  are  based 
on  the  frequency  dependence  of  the  observed  diffusion  length  and 
measurement  of  the  diffusion  velocity.  Such  measurements  yield 
both  the  diffusion  length  (L)  and  the  lifetime  (t),  which  are  related 
by  the  equation 


^^J.R.  Haynes  and  W.  Shockley,  Phys.Rev.  835  (1951). 
^“^A.A.  El-Difrawi,  T.H.  Yeh,  and  J.E.  Christopher,  J. 
Appl.  Phys.  £7,  272  (1976). 
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where  D is  the  diffusivity . The  mobility  (n)  and  diffusivity  are  in 
turn  related  by  the  Einstein  relation 

eD  = jxkT.  (2) 

Approximate  agreement  among  these  various  methods  was  found, 
and  on  this  basis  an  experimental  procedure  is  given  here  for  the 
measurement  of  both  the  bulk  diffusion  length  and  the  drift  mobility. 

6.2  EXPERIMENTAL  PROCEDURE 

The  measurement  of  drift  mobility  was  performed  using  a 
Jeol  JSM-Z  scanning  electron  microscope  modified  to  permit 
intensity -modulation  of  the  beam  by  deflection  past  an  aperture. 
Signal  recovery  of  the  phase  and  amplitude  information  was  by 
means  of  a network  analyzer.  The  experimental  design  is  illus- 
trated schematically  in  Figure  64.  A signal  generator  (HP8601) 
furnished  the  deflection  signal  as  well  as  a phase  and  amplitude 
reference  for  the  network  analyzer  (HP8407),  permitting 
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Figure  64,  Schematic  diagram  of  apparatus  for  narrowband  signal 
recovery  of  the  modulated  beam  response  in  amplitude  and 
phase . 
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measurements  in  the  range  of  100  kHz  to  110  MHz.  The  deflection 
’ signal  is  amplified  by  a 0.5W  broadband  amplifier  and  applied  to 

a resonant  deflection  structure.  A dc  bias  is  also  applied  to  the 
deflection  plates  to  yield  50%  duty  cycle.  The  deflection  plates  are 
0.5  X 2-in.  microstrip  substrates  placed  in  the  drift  space  beneath 
the  anode.  They  were  placed  as  far  apart  as  possible  to  minimize 
interference  with  the  beam.  Deflection  to  cut  off  a 25-kV  beam 
therefore  requires  some  200  volts.  The  plates  are  permanently 
left  in  place  and  do  not  interfere  with  conventional  operation  of  the 
SEM. 

The  sample  holder  was  provided  with  a MOSFET  preamplifier 
of  46  dB  gain  and  the  sample  stage  furnished  with  minature  50 -f} 
coax  to  the  chamber  wall.  Further  amplification  is  provided  out- 
side of  the  SEM.  To  increase  the  dynamic  range  beyond  the  limit 
imposed  by  the  10-kHz  noise  bandwidth  of  the  network  analyzer,  a 
[ crystal  filter  was  interposed  to  yield  a 10-Hz  noise  bandwidth. 

1 This  introduced  a phase  slope  of  90°  between  3-dB  points,  neces- 

sitating the  frequency  stabilization  of  the  signal  generator,  which 
was  characterized  by  100  Hz  of  incidental  frequency  modulation. 

A crystal-stabilized  oscillator  was  used  to  phase -lock  the  signal 
generator.  Resultant  phase  noise  was  then  on  the  order  of  0.3°. 
t (Network  analyzers  are  now  available  with  receiver  bandwidths  of 

10  Hz,  which  obviates  this  complication.) 

With  the  sensitivity  afforded  by  the  crystal  filter,  it  was  found 
that  residual  pickup  limited  the  dynamic  range.  In  excess  of  200  dB 
of  isolation  was  required  between  the  drive  and  output  circuitry. 

The  residual  pickup  was  simply  cancelled  b/  means  of  a signal  de- 
rived from  the  generator,  suitably  phase -shifted  and  attenuated. 

The  phase  and  amplitude  were  iteratively  adjusted  to  minimize  the 
' amplitude  of  the  observed  signal  and  to  obtain  a random  phase 

characteristic  of  pure  noise,  while  the  beam  current  was  turned  off. 

. Application  of  the  drift  field  to  the  sample  was  by  means  of  a 

pulse  generator  and  inverter,  which  could  be  adjusted  to  yield  a 
balanced  waveform  at  the  junction  position.  A small  (several  kT/e) 
. reverse  bias  was  applied  to  the  junction  to  reduce  the  junction 

capacitance  and  to  make  it  a more  ideal  "sink"  for  minority  car- 
riers, insensitive  to  small  bias  transients.  A typical  sample  is 
illustrated  in  Figure  65.  The  3 fl-cm  phosphorus -doped  chip  is 
furnished  with  two  large  electrodes  for  application  of  the  field,  and 
with  six  diffused  junctions  (of  1 -pm  depth)  in  the  central  area.  Two 
of  the  junctions  are  1-mil  x 20-mil  (under  the  long  metallization 
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Figure  65.  Silicon  chip  mounted  on  LSI  package  for  mobility 

measurement. 

fingers),  and  four  are  1-mil  x 1-mil  (at  the  ends  of  the  inter- 
mediate-length fingers).  Seven  ohmic  contacts  are  provided  at  the 
ends  of  the  small  fingers  for  potential  profiling.  One  of  the  bond- 
ing pads  is  furnished  with  small  fingers  for  calibration  of  the  mag- 
nification factor  . The  chip  is  mounted  on  the  insulated  substrate 
of  an  integrated  circuit  package. 

6.3  EXPERIMENTAL  RESULTS 

Preliminary  to  the  determination  of  mobility,  the  diffusion 
length  was  measured  under  ac -excitation  conditions  by  recording 
the  amplitude  response  as  the  beam  was  swept  in  a direction  nor- 
mal to  the  junction  periphery.  Modulation  frequencies  ranged 
from  30  kHz  to  10  MHz.  (The  30-kHz  data  were  taken  with  an 
Ithaco  353  logarithmic  phase -sensitive  detector.)  Data  obtained 
at  1 MHz  are  illustrated  in  Figure  66  where  results  are  shown  for 
a sweep  across  one  of  the  junctions.  The  phase  response  is  also 
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given.  A 25-keV  beam  was  used.  This  energy  represented  a com- 
promise between  considerations  of  ease  of  deflection  and  depth  of 
penetration,  which  governs  the  imporance  of  surface  effects. 

Beam  current  was  held  to  less  than  30  pA  to  minimize  the  volume 
over  which  high-injection  conditions  prevailed  (i.e.,  where  the 
excess  carrier  concentration  was  sufficient  to  modulate  the  local 
conductivity).  Since  the  pair  creation  energy  in  silicon  is  approxi- 
mately 3.6  5 eV,  each  25-keV  electron  generates  on  the  order  of 
7000  carrier  pairs,  so  that  high -injection  conditions  certainly 
prevail  within  the  energy  deposition  volume. 

An  analytical  expression  for  the  phase  and  amplitude  response 
for  excitation  at  angular  frequency  uj  has  been  given  by  McKelvey 
for  the  one -dimensional  case  :65 

V(x,uj)  = V(0,u))exp|-j^jexp|-ij^Y(UJ)).  (3) 

where  x is  the  distance  of  the  point  of  injection  from  the  junction. 
Here  L(uo)  is  the  ac  diffusion  length,  related  to  the  dc  diffusion 
length  L(0)  by  the  equation 


L(0) 


r,  / 

) 

1 + 1 1 + (u)  t)  1 

( 

( 

and  is  the  phase  factor  given  by 


(4) 


(5) 


The  amplitude  function  thus  yields  the  ac  diffusion  length  directly. 
The  data  of  Figure  66  exhibit  some  variation  in  the  apparent  local 
diffusion  length  and  this  is  due  to  deviations  from  one -dimensional 
behavior  near  the  junction,  as  well  as  to  the  effects  of  surface  re- 
combination. With  increasing  distance  from  the  junction,  the  dif- 
fusion length  approaches  the  bulk  value.  The  effects  of  surface 
recombination  have  been  treated  analytically  by  vanRoosbroeck, 
and  are  calculable  in  principle.  However,  explicit  fitting  of  the 

65 

J.P.  McKelvey,  Solid  State  and  Semiconductor  Physics 
(Harper  and  Row,  1966,  New  York),  pp.  439-440. 

66w.  vanRoosbroeck,  J.  Appl.  Phys . 380  (1955). 
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amplitude  function,  taking  into  account  surface  effects  and  detec- 
tor geometry,  is  deemed  to  be  intractable.  vanRoosbroeck  has 
shown  that  at  distances  from  the  junction  sufficiently  great,  sur- 
face recombination  affects  the  amplitude  only  as  a multiplicative 
constant.  It  is  therefore  preferable,  when  possible,  to  take  the 
experimental  approach  of  providing  sufficient  dynamic  range  to 
yield  the  bulk  value  of  diffusion  length  directly.  A first-order 
correction  for  deviation  from  one -dimensional  geometry  was 
made  (discussed  below)  and  a value  of  10.9  pm  for  the  diffusion 
length  obtained.  The  noise  level  of  Figure  66  corresponds  to 
about  4 dB  above  thermal  noise.  The  phase  characteristic  has  a 
range  of  150°  before  noise  dominates  and  indicates  a uniform 
signal  propagation  velocity.  Measurements  of  the  kind  shown  in 
Figure  66  have  previously  been  made  on  germanium  by  Munakata,^"^ 
at  lower  frequencies  and  for  longer  diffusion  lengths,  in  order  to 
verify  the  applicability  of  the  one -dimensional  model.  Consist- 
ency within  experimental  error  was  observed. 

The  drift  mobility  is  derived  most  directly  from  the  change 
in  signal  propagation  velocity  v,  as  observed  from  the  change  in 
phase  (0)  under  application  of  an  electric  field  (E).  We  have  then 


^drift  _ J. , . . _ UJx  / 1 1 \ 

E " E^'^E  ■ " e\0^(x)  ■ 0q(x)/ 


(6) 


Vq  is  governed  strictly  by  diffusion  whereas  v^  has  the  drift  veloc- 
ity superimposed  on  it.  Field  strengths  were  as  high  as  30  to 
50  V /cm  and  short  pulses  were  employed  to  minimize  sample 
heating.  Use  of  the  crystal  filter,  however,  limited  the  signal 
bandwidth  to  such  an  extent  that  the  phase  did  not  settle  within  the 
desired  pulse  lengths.  Equivalent  information  was  therefore 
derived  from  the  amplitude  function.  Figure  67  shows  the  ampli- 
tude profiles  under  zero  field  as  well  as  under  retarding  and  ac- 
celerating fields.  Analysis  assumes  that  equal  amplitude  signals 
correspond  to  equal  propagation  times.  Thus, 


(7) 


C.  Munakata  and  T.E.  Everhart,  Japan  J.  Appl.  Phys . Ml 
913  (1972). 
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Figure  67.  The  dependence  of  amplitude  on  distance  from  the  junc- 
tion under  field-free  (A,  □ ) conditions,  as  well  as  accelerating 
pulsed  field  (O),  and  retarding  pulsed  field  (0)  conditions,  for  a 
1 -MHz  modulated  beam.  The  applied  field  was  11  V/cm.  The 
corresponding  sets  of  curves  have  been  separated  vertically  for 
clarity,  so  the  amplitude  scale  is  arbitrary. 

where  xg  and  Xq  refer  to  equal -amplitude  coordinates  on  the  respec- 
tive curves.  Mobilities  determined  in  this  maxuier  were  250  cm^/V-sec 
and  270  cm^/V-sec  on  the  two  sides  of  the  junction.  The  curves 
shown  in  the  figure  should  intersect  at  the  junction  periphery  and 
the  deviation  may  be  due  to  conductivity  modulation  near  the  point 
of  injection. 

Mobilities  were  also  measured  using  phase  information,  with- 
out the  crystal  filter.  Typical  values  were  consistent  with  those 
found  above,  centering  about  240  cm^/V-sec.  These  values  are  low 
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compared  to  the  mobility  of  about  450  cm^/V-sec  expected  for  holes 
in  3 fJ-cm  n-type  material,  (Possible  reasons  for  the  discrepancy 
will  be  discussed  below.  ) The  finding  of  an  unexpected  low  mo- 
bility motivated  a comparison  with  other,  field -free,  methods  of 
determining  the  mobility. 

The  mobility  may  be  determined  from  the  diffusivity  using 
the  Einstein  relationship  [Eq.  (Z)  ] if  both  the  lifetime  and  diffusion 
length  are  known.  Both  of  these  quantities  may  be  determined 
from  measurement  of  L(u!  ) or  the  signal  propagation  velocity  at  a 
number  of  frequencies,  encompassing  such  frequencies  that 
U)  T » 1.  (See  Eqs . (4)  and  (5).)  We  have  obtained  a mobility  value 
by  comparison  of  the  1-MHz  diffusion  length  with  the  dc  value,  as 
well  as  by  measurement  of  the  signal  propagation  velocity  for  a 
wide  range  of  frequencies. 

Results  for  the  measurement  of  the  dc  diffusion  length  are 
shown  in  Figure  68.  The  experimental  method  for  this  measure- 
ment has  been  described  elsewhere.^®  Since  the  noise  limit  of 
the  method  is  imposed  by  the  1/f-noise  of  the  (unpassivated)  detec- 
tor junction,  dynamic  range  equivalent  to  that  of  the  ac  measure- 
ment is  obtained  by  means  of  signal  averaging.  The  noise  level  in 
this  measurement  is  on  the  order  of  a few  pA.  The  amplitude  func- 
tion is  qualitatively  similar  to  that  of  Figure  66  in  that  the  apparent 
diffusion  length  is  foreshortened  near  the  junction,  due  to  surface 
effects  and  geometric  factors,  and  approaches  a constant  value  at 
some  distance  from  the  junction.  Since  the  apparent  diffusion  length 
is  spatially  varying,  an  improvement  on  a one -dimensional  treat- 
ment is  desirable.  The  steady-state  solution  of  the  diffusion  equa- 
tion in  spherical  coordinates,  for  point  source  excitation,  is  given 
by 

n(r)  = — exp(-r/L)  , (8) 

r 

where  n is  the  concentration  at  distance  r from  the  source  and  A 
is  an  assumed  source  strength.  The  measured  signal  is  directly 
proportional  to  the  value  of  n at  the  junction.  A fit  of  the  functional 
form  (1/x)  exp  (-x/L)  to  the  data  yielded  a spatially  invariant  value 
of  diffusion  length  of  17.6  pm,  which  was  taken  to  be  the  bulk  value. 
The  same  considerations  apply  to  the  1-MHz  data.  Verification 
that  the  data  for  the  two  cases  are  analogous  is  provided  in  a log -log 

68s.  Othmer  and  O.L.  Curtis,  Jr.,  IEEE  Trans.  Nucl.  Sci. 

204  (Dec.  1973). 
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Figure  68.  Measurement  of  the  minority-carrier  current  as  a 

function  of  distance  from  the  point  of  injection  to  the  jxinction, 
for  the  determination  of  the  dc  diffusion  length.  Shown  is  the 
output  of  a multichannel  scaler  after  summation  of  traces  for 
signal -to -noise  enhancement.  The  amplitude  scale  is 
arbitrary. 

plot  of  the  amplitude  function,  .shown  in  Figure  69.  The  data  sets 
have  been  superimposed  by  translation  along  the  two  coordinate 
axes,  and  the  correspondence  demonstrates  that  the  only  difference 
between  them  is  the  spatial  scale  (i.e.,  the  bulk  diffusion  length). 

The  same  fvuictional  form  was  therefore  used  to  fit  these  data  and 
the  value  of  10.9  pm  already  referred  to  was  obtained.  Using  Eq.  (4), 
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Figure  69.  Comparison  of  the  dc  diffusion  length  and  the  1-MHz 
diffusion  length.  The  scales  are  arbitrary  and  the  curves 
have  been  displaced  along  the  coordinate  axes  for  superposition. 


the  lifetime  is  determined  to  be  640  ns,  the  diffusivity  4.84  cm^/ 
sec,  and  the  mobility  190  cm^/V-sec.  The  precision  of  this  value 
is  poor,  due  to  the  use  of  the  rather  low  frequency  of  1 MHz  for 
this  short  lifetime  (out  = 4).  An  error  of  10%  in  the  1 -MHz  diffusion 
length  implies  an  error  of  some  26%  in  mobility  in  this  case.  Never- 
theless, the  low  values  of  mobility  found  above  with  application  of 
an  electric  field  stand  confirmed  to  within  the  accuracy  of  this 
de  te  r mination . 

The  second  field-free  method  of  determining  the  drift  velocity 
is  based  on  phase  measurements  of  the  signal  propagation  velocity. 
The  phase  shift  in  Eq.  (3)  has  the  form  yx/L(0),  so  that  the  signal 
propagation  velocity  is 

V = L(0)  (uj/y),  (9) 


from  which  the  lifetime  may  be  determined  if  the  dc  diffusion 
length  is  also  known.  The  signal  propagation  velocity  is  closely 
related  to  the  diffusion  velocity,  given  by  L(0)/t.  In  the  low- 
frequency  limit,  the  signal  propagation  is  equal  to  twice  the 
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diffusion  velocity.  The  data  obtained  for  signal  propagation  veloc- 
ity are  shown  in  Figure  70  for  all  of  the  measurement  frequencies. 
Only  the  1-MHz  and  3. 6 -MHz  data  points  were  obtained  with  a 
crystal  filter.  Precision  of  the  other  data  points  is  considerably 
poorer  - at  high  frequencies  due  to  limited  dynamic  range,  and  at 
low  frequencies  due  to  the  fact  that  the  total  phase  change  measured 
was  rather  small.  The  curve  shown  is  a fit  to  the  data  which 
weights  the  two  best  data  points  most  heavily.  This  yields  a life- 
time of  320  ns,  and  thus  a diffusivity  of  9.7  cm^/sec  or  a mobility 
of  370  cm^/V-  sec.  Since  the  determination  is  based  on  data  taken 


at  such  low  frequencies,  the  margin  of  error  could  easily  be  +100/ 
-70  cm^/Y-  sec  in  this  determination.  However,  the  general  con- 


formance of  the  data  to  expectations  tends  to  justify  the  use  of  a 


one -dimensional  model,  and  the  use  of  this  method  of  determining 


the  drift  mobility. 


f (MHz) 


Figure  70.  Plot  of  signal  propagation  velocity,  determined  from 
phase  measurements,  vs  modulation  frequency.  The  fit 
shown  is  based  on  the  one -dimensional  model.  Refer  to 
the  text  for  discussion  of  the  fit. 
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6.4  DISCUSSION 


Comparison  has  been  made  between  methods  of  determining 
drift  mobility  which  do  not  require  application  of  a field  and  those 
using  a pulsed  electric  field.  The  data  are  recapitulated  in 
Table  5,  Considerable  spread  in  the  measured  values  was  found 
among  the  various  techniques,  but  a general  tendency  for  the  values 
to  be  low  compared  to  expectations  was  manifested.  This  is  sur- 
prising, and  several  possible  explanations  have  been  considered. 
First,  under  high-injection  conditions,  transport  is  governed  not 
by  the  drift  mobility,  but  rather  by  the  ambipolar  mobility,  which 
tends  to  zero  in  the  intrinsic  limit.  At  the  low  beam  currents  em- 
ployed here,  such  high-injection  conditions  were  determined  to 
exist  only  within  the  energy  deposition  volume.  Therefore,  it  is 
unlikely  that  this  factor  is  principally  responsible  for  the  observa- 
tion of  a low  mobility.  Moreover,  the  field -free  methods  are 
governed  under  high -injection  conditions  by  the  ambipolar  diffu- 
sivity,  which  does  not  vanish  in  the  intrinsic  limit,  and  is,  in  fact, 
larger  than  the  low-injection  hole  diffusivity  in  the  case  of  n-type 
samples.  A second  possible  explanation  is  that  the  mobility  is  re- 
duced due  to  local  heating,  but  this  factor  has  likewise  been  found 
insufficient  to  explain  the  observed  discrepancy.  Finally,  it  is 
conceivable  that  the  present  measurements  are  inordinately  influ- 
enced by  the  presence  of  the  surface.  An  upper  limit  to  the  pos- 
sible influence  of  the  surface  is  indicated  by  the  channel  mobilities 

Table  5. 

Mobilities  Determined  by  Various  Methods 


Measurement 

Mobility 

Technique 

(cm  /V-sec) 

Pulsed  Field, 

Amplitude  Function 

250-270 

Pulsed  Field, 

Phase  Function 

240 

ac -Diffusion  Length 

190 

Signal  Velocity 

370 

(Expected  Value) 

(450) 
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of  MOS  transistors,  where  transport  is  influenced  by  the  interface 
and  in  which  mobilities  on  the  order  of  one -half  of  bulk  mobilities 
are  observed.  This  effect  would  seem  to  be  insufficient  to  explain 
the  discrepancy  in  the  present  case,  which  approaches  the  condi- 
tions for  a bulk  measurement  much  more  closely. 

The  spread  in  the  observed  mobility  values  among  the  various 
techniques  is  ascribed  principally  to  the  nonoptimal  choice  of  oper- 
ating frequencies  in  the  case  of  the  field -free  techniques,  and  to 
the  fact  that  mobility  determinations  become  vastly  more  difficult 
at  short  diffusion  lengths.  Radiation  effects  studies  at  such  short 
diffusion  lengths  therefore  appear  to  be  unattractive.  On  the  other 
hand,  the  prospect  of  using  these  methods  for  the  characterization 
of  solar  cells  appears  promising.  Silicon  solar  cells  must  exhibit 
diffusion  lengths  of  at  least  50  (j.m  in  order  to  have  any  prospect  of 
economic  viability  and  the  above  techniques  should  be  capable  of 
yielding  good  mobility  values  in  the  range  of  diffusion  lengths  above 
20  pm. 

The  principal  shortcomings  of  the  above  techniques  relate  to 
the  influence  of  the  surface  on  the  measurements  and  the  assump- 
tion of  one -dimensional  character  of  transport.  Both  of  these  lia- 
bilities are  remedied  if  measurements  of  the  ac  diffusion  length 
are  made  with  a bevel  geometry,  in  which  the  bombarded  surface 
and  the  collecting  surface  form  a wedge  of  arbitrary  apex  angle 
(on  the  order  of  10°).  Transport  then  occurs  through  the  bulk  and 
the  effect  of  the  surface  on  the  measurements  is  totally  negligible 
as  long  as  the  surface  recombination  velocity  at  the  injecting  sur- 
face is  uniform.  The  transport  is  also  more  one -dimensional  in 
character  by  virtue  of  the  fact  that  the  diffusion  front  is  tangential 
to  the  collecting  surface.  DC  measurements  of  diffusion  length  by 
such  methods  have  been  quite  common,  particularly  for  GaAs . 

The  simple  expedient  of  measuring  the  amplitude  and  phase  response 
using  two  appropriate  frequencies  should  yield  both  the  dc  diffusion 
length  and  the  drift  mobility  for  the  material  as  well  as  a measure 
of  the  spatial  uniformity  of  these  parameters.  One  frequency  should 
be  chosen  such  that  ujt  « 1,  to  yield  a diffusion  length  close  to  the 
dc  value,  and  the  other  should  be  chosen  such  that  lUT  > 10,  to  give 
values  for  the  lifetime  on  the  basis  of  both  the  phase  and  the  ampli- 
tude response.  Jointly,  then,  these  parameters  yield  the  mobility. 
The  determination  of  good  mobility  values  by  this  method  requires 

^9w.H.  Hackett,  Jr.,  J.  Appl.  Phys . 43,  1649  (1972). 
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precision  not  usually  associated  with  diffusion  length  measure- 
ments. We  have  used  the  bevel  geometry  for  determining  diffusion 
lengths  of  40  pm  on  solar  cells  with  5%  precision,  and  have  found 
agreement  of  this  value  with  a short-circuit  current  measurement 
of  the  diffusion  length  using  the  uniform  generation  of  a Cobalt-60 
gamma  source.  Measurement  with  the  bevel  geometry  clearly  does 
afford  the  desired  accuracy. 


^ t- 
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Washington,  DC  20362 

Andrews  AFB 

Attn;  SEA-9931,  S.  Barham 

Washington,  DC  20334 

5EA-04531 

Atm:  DLCA 

Officer  -in -Charge 

Air  Force  Weapons  Laboratory 

Naval  Surface  Weapons  Center 

KirUand  AFB.  NM  87117 

WTiite  Oak,  Silver  Spring,  MD  20910 

Atm:  ELP,  TREE  Section 

Atm:  Code  F3  1 

ELA 

Code  F30 

ELXT 

EL,  C.  Baum 

Commander 

ELS 

Naval  Surface  Weapons  Center 

NXS 

Dahlgren  Laboratory 

J.  Mullis 

Dahlgren.VA  22448 

NT 

Attn:  DG-50,  W.  Holt 

SUL 

DEX 

Commander 

Naval  Weapons  Center 

R.  Maier 

China  Lake,  CA  93555 

Headquarters 

Attn;  Code  533,  Tech.  Lib. 

Electronic  Systems  Division  IN 

Department  of  the  Air  Force 

Commanding  Officer 

HanscomAFB,  MA  01731 

Nuclear  Weapons  Tng  Center  Pacific 
Department  of  the  Navy 

Atm:  rNDC/21 

Naval  Air  Station,  North  Island 

Headquarters 

San  Diego,  CA  92135 

Electronic  Systems  Division  ■ YSEA 

Atm:  Code  32 

Department  of  the  Air  Force 

Hanscom  AFB,  MA  01731 

Office  of  Naval  Research 

Arlington,  VA  22217 

Atm;  YSEA 

Attn:  Code  427 

Code  421,  Doran  W.  Padgett 
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Assistant  to  the  Secretary  ot 
Defense 
Atomic  Energy 
Washington.  DC  20301 
Attn:  Executive  Assistant 

Director 

Command  ic  Control  Tech  Center 
Department  of  Defense 
The  Pentagon,  Rm  BE  633 
Washington,  DC  20301 
Attn:  C-362,  Mr.  Adkins 


Defense  Communication  Engineer 
Center 

1360  Wiehle  Avenue 
Reston,  VA  22090 
Attn:  Code  R410 

Defense  Documentation  Center 
Cameron  Station 
.Alexandria,  V’A  22314 
Attn:  DD  (12) 

DDC-TCA  (2) 

Ccrrunander 

Defense  Electronic  Supply  Center 
1307  Wilmington  Pike 
Dayton.  OH  43401 

(Thru  DSA  for  CXWDI) 

Attn:  D ESC -ECS 

Dir  ector 

Defense  Intelligence  Agency 
Washington.  DC  20301 
Attn:  DB-4C 


Director 
Defense  Nuclear  Agency 

I Washington,  DC  20305 

[ Attn:  STRA 

I STVL 

I RAEV 

[ DOST 

I TITL  (4) 

[Commander 

Field  Command 
Defense  Nuclear  Agency 
Kirdand  AFB,  NM  37115 
Attn:  FCPR 

Director 

Interser*/tce  Nuclear  W'eapons 
School 

Kirtiand  AFB.  NM  87115 
Attn:  TTV 

Director 

Joint  Strat  TGT  Planning  Staff 
Offutt  AFB 
Omaha,  68113 
Attn:  JLTW-2 


Chief 

Livermore  Division  Fid  Command 
DNA 

Lawrence  Livermore  Laboratory 
P.O.  Box  308 
Livermore,  CA  94550 
Attn:  FCPRL 

National  Communications  System 
Office  of  the  Manager 
Department  of  Defense 
Washington,  DC  20305 
Attn:  NCS-TS,  C.  Bodson 

Director 

National  Security  Agency 
Department  of  Defense 
Fort  George  G,  Meade,  MD  20755 
AtGi:  R-52,  O.  Van  Cunten 

Under  Secy  of  Def  for  Research 
Engineering 
Department  of  Defense 
Washington,  DC  20301 
Attn:  Strategic  i Space  Systems  (OS) 

Director 

3MD  Advanced  Technology  Center 
Hunts%'ille  Office 
Department  of  the  Arm.y 
P.O.  Box  1500 
Huntsville,  AL  35307 
Attn:  ATC-T 
F.  Hoke 

Commander 

BMD  Systems  Command 
Department  of  the  .Army 
P.O.  Box  1500 
Huntsville,  AL  35807 
Attn:  BMDSC-TEN,  N.  Hurst 
BMDSC-HEN 

Commander 

ERADCOM  Technical  Support  Activit> 
Department  of  the  Army 
Fort  Monmouth,  NJ  07703 
Attn:  DRSEL-TL-EN 

DRSEL-PL-ENV,  H.  Bomke 
A.  Cohen 

DR5EL-TL-MD.  G,  Gaule 
DRSEL-NL-RC,  R,  BrOv'.Ti 

Commander 
Fort  Huachuca 
Department  of  the  Army 
Fort  Huachuca,  AE  35613 
Attn;  Tech  Ref  Div 

Commander 

Redstone  Scientific  Info  Center 
U.S.  .Army  RitD  Command 
Redstone  Arsenal,  AL  35809 
•Atto:  Chief,  Documents  (3) 
DRDMI-TBD 


Commander 

Harry  Diamond  Laboratories 
Department  of  the  Army 
2800  Powder  Mill  Road 
Adelphi,  MD  20783 

Attn:  J.  Thompkins 

DELKD-NP,  F.  Wimenitz 
DELHD-NP 

DRXDO-EM-1,  J.  .Miletta 
DELHD-EM,  R.  .McCoskey 
DELHD-RBA,  J.  Rosado 
DELHD-RBG 

DELHD-RBH,  F.  .McLean  (10) 
DELHD-TF,  R.  Oswald,  Jr. 
DELHD-TI  (Tech  Lib)  (3) 
DELHD-RB,J.  McGarrity 
DELHD-TI  ^Tech  Reports) 
Chairman,  Editorial  Comm. 
DRXDO-AS  (Records) 


Secretary  of  the  Army 
Washington,  DC  20301 
Attn;  ODUSA  or  Daniel  Willard 

Commander 

U.S.  Army  Armament  Research  i: 

Development  Command 
Dover,  NJ  07801 
•Attn;  DRDAR-LCN-F 
DRDAR-LCA-PD 
DRDAR-TSS  No.  5^ 

Director 

U.S.  Army  Ballistic  Research  Labs 
Aberdeen  Proving  Ground,  MD  21005 
Atm:  DROAR-ELT 

DRXBR-AM,  W.  VananfA’erp 
DRXRD-BVL.  D.  Rigotti 
STEAP-RF,  R.  Harrison 

Chief 

U.S.  .Army  Communications  Sys  Agcy 
Fort  Monmouth,  7'J  07793 
Attn:  SCCM-AD-SV  (Library) 
CCM-RD-T,  S.Krevsky 

Commander -in -Chief 
U.S.  .Army  Europe  Seve.nth  .Army 
APO  New  York  09403 
I'Heidelberg ) 

Atm;  CDSCE-E  AEAGE-PI 
Commander 

U.S.  Army  Intelligence  & Sec  Cmd 
Arlington  Hall  Station 
4000  Arlington  Blvd. 

Arlington,  VA  22212 

.Atm:  lARDA-OS,  R.  Burkhardt 
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Headquarters 

Electronic  Systems  Division,  AFSC 
Hanscom  AFB.  MA  01731 
Attn:  Technical  Library 

Commander 

Foreign  Technology  Di'/ision,  AFSC 
Wright-Patterson  AFB,  OH  45433 
Attn:  ETDP 

ETDP,  B.  Ballard 

Commander 

Rome  Air  Development  Center, 
Griffisa  AFB,  NY  13441  AFSC 

Attn;  RADC/RBRP,  C.  X-ane 

RADC/RBRAC,  I.  L.  Krulac 

Commander 

Rome  Air  Development  Center, 
Hanscom  AFB,  MA  01731  AFSC 

Attn:  £SR,  B.  Buchanan 
ES£,  A.  Kahan 
ETS,  R.  Dolan 

Space  Jit  Missile  Systems  Organiza- 
tion/DY 

Air  Force  Systems  Command 
P.  O.  Box  92960 
Worldway  Postal  Center 
Los  Angeles,  CA  90009 
(Technology) 

Attn:  DYS 

Space  &{  Missile  Systems  Organiza- 
tion /IN 

Air  Force  Systems  Command 
P.  O.  Box  52960 
Worldway  Postal  Center 
Los  Anaeles,  CA  90009 
(Intelligence) 

Attn:  IND 


Space  Missile  Systems  Organiza- 

tion/SK 

Air  Force  Systems  Command 
P.  O.  Box  92960 
Worldway  Postal  Center 
Los  Angeles,  CA  90009 
(Space  Comm  Systems) 

Attn:  SKP.  P.  Sudler 

Space  Sc  Missile  Systems  Organiza* 
tion/YA 

Air  Force  Systems  Contmand 
P.  O.  Box  92960 
Worldway  Postal  Center 
Los  Angeles,  CA  90009 
Attn:  YAS 

Director  (INWS) 

34i6th  Technical  Training  Squadron 
Air  Training  Command 
Department  of  the  Air  Force 
Kirtland  AFB,  NM  87115 
Attn:  TTV 

Department  of  Energy 
Albuquerque  Operations  Office 
P.  O.  Box  5400 
Albuquerque,  NM  87115 
Attn:  Doc  Con  for  WSSB/OSD, 

R.  Shay 

University  of  California 
Lawrence  Livermore  Laboratory 
P.  O.  Box  808 
Livermore,  CA  94550 
Atm:  Doc  Con  for: 

L-156,  R.  Kalibjian 
L-389,  R.  Ott 

L-96,  H.  Kruger  (Class  L-'- 
L-125,  J.  KeUer 
Technical  Information  Dept. 


Aerojet  Electro -Systems  Co. 

Div.  of  Aerojet-General  Corp. 

P.  O.  Box  296 
Azusa.  CA  91702 
Attn:  SV/87n/70 

Aerospace  Corp. 

P.  O.  Box  92957 
Los  Angeles,  CA  90009 
Attn:  V.  Josephson 

I.  Garfunkel 
W.  Willis 
R.  Crolius 

J.  Reinheimer 

AVCO  Research  it  Systems  Group 
201  Lowell  Street 
Wilmington,  MA  01887 
Attn:  W,  Eroding 

Battelle  Memorial  Institute 
505  King  Avenue 
Columbus,  OH  43201 
Attn:  R.  Thatcher 
R.  Blazek 
D.  Hamman 

BOM  Corp. 

P.  O.  Box  9274 
Albuquerque  International 
Albuquerque,  NM  37119 
Attn:  D.  Alexander 
Marketing 

Bendix  Corp. 

Communication  Division 
East  Joppa  Road 
Baltimore,  MD  212C-* 

Attn:  Document  Control 

Bendix  Corp. 


Space  i:  Missile  Systems  Organiza- 
tion /MN 

Air  Force  Systems  Command 
Norton  AFB.  CA  92409 
(Minuteman) 

Attn:  MNNH 
MNNG 

Space  U Missile  Systems  Organ!  za- 
tion/RS 

Air  Force  Systems  Command 
P.  O.  Box  92960 
Worldway  Postal  Center 
Los  Angeles,  CA  90009 
(Reentry  Systems) 

Attn:  RSSE,  LTC  K.  Gilbert 
RSMA,  Lt.  D.  Higgins 
RSMC,  Capt.  Collier 


Office  of  Military  Application 
Department  of  Energy 
Washington,  DC  20545 
Attn;  Doc  Con  for  Classified  Library 


Research  Laboratories  Di\'ision 
Bendix  Center 
Southfield,  MI  48075 
Attn:  M,  Frank 


Sandia  Laboratories 
P.  O.  Box  5800 
Albuquerque,  NM  87115 
Atm:  Doc  Con  for: 

3141 
J.  Hood 
F • Coppage 
B.  Gregory 


Boeing  Co. 

P.  O.'box  3707 
SeatUe,  WA  9 3124 
Atm:  8K-38 

Aerospace  Library 
I.  Arimura 
R.  Caldwell 
C.  Rosenberg 


Central  Intelligence  Agency 

Washington,  DC  20505 

Atm;  RD’Sl,  Rm  3G48,  HQ  Bldg. 


Boo 2 -Allen  and  Hamilton,  Lnc. 
776  Shrewsbury  Avenue 
Tinton  Falls,  NJ  07724 
Atm:  R.  Chrisner 


Strategic  Air  Command /XPFS 
Offutt  AFB.  , NB  68113 
Attn:  XPFS,  M.  Carra 

NRI-STINFO  Library 


Department  of  Transportation 
Federal  Aviation  Administration 
Headquarters  Sec  Div.,  ASE-300 
800  Independence  Ave.,  SW 
Washington,  DC  20591 
Atm:  ARD-350 


Brown  Engineering  Company,  Lnc. 
Cummings  Research  Park 
Huntsville,  AL  35807 
Atm:  J.  MeSwain 
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Cdlifornia  Institute  of  Technology 
Jet  Propulsion  Lab 
4800  Oak  Grove  Drive 
Pasadena,  CA  91103 
Attn:  J.  Bryden 
A.  Stanley 

Charles  Stark  Draper  Lab.,  Inc. 
555  Tech.nology  Square 
Cambridge,  MA  02139 
.A.ttn:  K.  Ferrig 
P.  Kelly 

R.  Halrmaier 
P.Greiff 
H.  Hueramler 

Computer  Sciences  Corp. 

1400  San  .Mateo  Blvd..  5F 
.Albuquerque,  NM  87108 
.\ttn:  A.Schiif 

Control  Data  Corp 
P.  O. Box  C 

Mi.nneapolis,  MN  5 5440 
.Attn:  J.  Meehan 

Cutler  •Ha.mmer  , Lnc  . 

.AIL  Di’/ision 

Comae  Road 

Deer  Park,  .\Y  11729 

.\ttn:  Cent  Tech  Files,  .A.  Anthony 

Dikewood  Lndustries,  Inc. 

1009  Bradbury  Drive,  SE 
.Albuequerque,  7.'M  87106 
.Attn:  L . Davis 

E-Systen'.s,  Inc. 

ECI  Di’nsion 
P.  O.  Box  12248. 

St.  Petersburg,  FL  33733 
Attn:  R . F re.nch 

F-Systems,  Inc. 

Greenville  Division 
P . O . Box  1 03o 
Greenville,  TX  7 5401 
.•\,ttn:  Division  Library 

Effects  Technology,  Inc. 

5333  Hollister  .Avenue 
Santa  Barbara,  C.A  931i] 

Attn:  E.  Steele 

EX -Cal,  Lnc . 

First  National  Bldg.,  E 
Suite  1316 

.Albuquerque,  \M  87108 
.Attn:  .R.  Dickhaut 


Fairchild  Industries,  Inc, 

Sherman  Fairchild  Tech,  Center 
20301  Century  Blvd. 

Germantown,  .MD  20767 

Attn:  Mgr  Config  Data  &c  Standards 

University  of  Florida 
Attn:  Security  Officer 
P.  O.  Box  284 
Gainesville,  FL  32601 
Attn;  H.  Sisler 

Ford  Aerospace  Communications 
Ford  &£  Jamboree  Roads 
Newport  Beach,  CA  92663 
Attn;  E.  Poncelet,  Jr. 

K.  Attinger 

Tech  Info  Services 

Ford  Aerospace  ic  Communications 
3939  Fabian  Way 
Palo  Alto,  CA  94303 
Attn:  Technical  Library 

S.  Crawford 

General  Electric  Co. 

Space  Division 
Valley  Forge  Space  Center 
P . O . Box  85  5 5 
Philadelphia,  PA  19101 
.Attn:  D.  Tasca 
D.  Long 

L.  Jeffers 
L.  Chasen 

J.  Peden  VFSC,  Rm.  4230M 
L.  Sivo 
J.  .Andrews 

General  Electric  Co. 

Re-Entry  k Environmental  Systems 
P.  O.  Box  7722 
Philadelphia,  PA  19101 
Attn:  W.  Patterson 

J.  Palchefsky,  Jr. 

Technical  Library 
R.  Benedict 

General  Electric  Co. 

Ordnance  Systems 
100  Plastics  Avenue 
Pittsfield.  MA  01201 
Attn:  J.  Reidl 

General  Electric  Co. 

Aerospace  Electronics  Systems 
French  Road 
Utica.  NY  13503 
Attn:  C.  Hewison 


General  Electric  Co.  - Tempo  Cent. 
: airchild  Camera  k Instrument  Corp  for  Advanced  Studies 

464  Ellis  Street  P.  O.  Drawer  QQ 

Mo'untaLn  View,  C.A  94040  Santa  Barbara.  CA  93102 

Att.n:  Sec  Con  for  D.  Myers  .Attn-  R.  Rutherford 

.M.  Espig 
DAS  LAC 


General  Electric  Co.  - Tempo 
Alexandria  Office 
Huntington  Bldg.,  Suite  300 
2560  Huntington  Ave. 

Alexandra,  VA  22303 
Attn;  DASIAC 

General  Research  Corp. 

Santa  Barbara  Division 
P.  O.  Box  6770 
Santa  Barbara,  CA  93111 
Attn:  Tech  Information  Office 

Georgia  Institute  of  Technology 
Georgia  Tech  Research  Lnscitute 
Atlanta,  CA  30332 
Attn;  R.  Curry 

Georgia  Institute  of  Technology 
Office  of  Contract  Administration 
•Attn:  Rsch  Security  Coordinator 
.Atlanta,  GA  30332 
Attn;  Res  it  Sec  Coord  for  H.  De.nr.y 

Grumman  -Aerospace  Corp. 

S.  Oyster  Bay  Road 
Bethpage,  NY  11714 
Attn;  J.  Rogers 

GTE  Sylvar.ia,  Inc.  - Electronics  Sys. 

Group  - Eastern  Division 
77  A Street 
Needham,  MA  02194 
Attn:  L.  Blaisdell 

C.  Thornhill,  Librarian 

GTE  Sylvania,  Inc. 

189  B Street 

Needham  Heights , MA  02194 
Attn:  H&V  Group 
J.  Waldron 
P.  Fredrickson 

Harris  Corp. 

Electronics  Systems  Division 
P.  0.  Box  37 
Melbourne,  FL  32901 
Attn:  C.  Davis 

Harris  Corp. 

Harris  Semiconductor  DiNosion 
P.  O.  Box  833 
Melbourne.  FL  32901 
Attn:  Mngr  Linear  Eng 

T .  Clark 

Mngr  Bipolar  Digital  Eng 

Hazeltine  Corp. 

Pulaski  Road 

Greenlawn,  NY  11740 

Attn:  Tech  Info  Ctr,  .M.  Waite 

Honeywell,  Inc. 

Avionics  Division 
2600  Ridgeway  Parkwav 
Minneapolis,  MN  55413 
Attn:  R.  Johnson 
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Honeywell,  Inc. 

Avionics  Division 
13350  U.S.  Highway  19,  N. 

St.  Petersburg,  FL  33733 
Attn:  MS  725-5 

Honeywell,  Inc. 

Radiation  Center 
2 Forbes  Road 
Lexington,  MA  02173 
Attn:  Technical  Library 

Hughes  Aircraft  Co. 

Centinela  and  Teale 
Culver  City,  CA  90230 
Attn:  K.  Walker 

CTDC,  o/EllO 
D.  Binder 

Hughes  Aircraft  Co. 

ITl  Segundo  Site 
P.  O.  Box  92919 
Los  Angeles,  CA  90009 
Attn;  E.  Smith,  MS  A620 

IBM  Corp. 

Route  17C 
Cwego,  NY  13827 

Attn:  Electromagnetic  Compatability 
Mono  Memory  Systems 

Institute  for  Defense  Analyses 
400  Army -Navy  Drive 
Arlington,  VA  22202 
Att.n:  Tech  Info  Ser-/ices 

Lon  Physics  Corp. 

S.  Bedfo rd  Street 
Burlington,  MA  01803 
Attn:  R.  Evans 

LRT  Corp. 

P.  0.  Box  81087 
San  Diego,  CA  92133 
.A,ttn:  MDC 

a.  Mertz 

Systems  Effects  Div. 

Physics  Div. 

JAYCOR 

205  S.  ^Vhiting  Street,  Suite  500 
Alexandria,  VA  22304 
Attn:  R.  Sullivan 

C.  Turesko 

Kaman  Sciences  Corp. 

P.  O.  Box  7463 
Colorado  Springs,  CO  80933 
Attn:  Dir  Science  4c  Tech  Div 
W,  Rich 
President 

D.  Bryce 
W.  Ware 
J.  Lubell 


Litton  Systems  , Lnc.  - Guidance 
and  Control  Systems  Div. 

5500  Canoga  Avenue 
WoodUndHUls,  CA  91364 
Attn:  VV.  .Mras 
J.  Retzler 

Lockheed  Missiles  & Space  Co.,  Inc 
P.  O.  Box  504 
Sunnyvale,  CA  94086 
Attn:  E.  Smith 

B.  Kimura 
L.  Rossi 

S.  Taimuty 
D.  Wolfhard 

Lockheed  Missiles  Space  Co.,  Inc 
3251  Hanover  Street 
Palo  .Aito,  CA  94304 
•Attn:  Tech  Info  Ctr,  D/Coll 
Reports  Library 

M.I.T.  Lincoln  Lab 
P.  O.  Box  73 
Lexington,  MA  02173 
Attn:  L.  Loughlin 

Library,  A-082 

Martin  Marietta  Corp. 

Orlando  Division 
P.  O.  Box  3837 
Orlando,  FL  32303 
Attn:  TlC/'MP-30 

Martin  Marietta  Corp. 

Denver  Division 
P.  O.  Box  179 
Denver,  CO  80201 
Attn:  P.  Kase 
Goodwin 

Research  Library 

•McDonnell  Douglas  Corp. 

P.  O.  Box  516 

St,  Louis,  NIC  63166 

Attn:  Library 

T.  Ender 

McDonnell  Douglas  Corp. 

5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 
Attn:  P.  .Albrecht 

Mission  Research  Corp. 

P . O . Drawer  719 
S<tnta  Barbara,  CA  93102 
Attn:  M.  Van  Blaricum 

Mission  Research  Corp. 

EM  System  Applications  Div. 

1400  San  Mateo  Blvd.  SE,  Suite  A 
Albuquerque,  NM  37  108 
Attn:  David  E.  Merewether 


Mission  Research  Corp.  - San  Diego 
P.  O.  Box  1209 
La  JoUa,  CA  92038 
Attn:  J.  Raymond 

V.  A.  J.  Van  Lint 

National  Academy  o:  Sciences 
, Committee  on  Atmospheric  Sciences 
2101  Constitution  Avenue,  NV/ 
Washington,  DC  20418 
Attn:  R. Shane 

University  of  New  .Mexico 
Electrical  Eng  h Computer  Science 
Department 

-Albuquerque,  NM  87131 
, Attn:  Harold  Southward 

Northrop  Corporation 
Northrop  Research  u Tech  Center 
One  Research  Park 
Palos  Verdes  Peninsula,  CA  90274 
Attn:  J.  R.  Srour 
S.  Othmer 
Library 

Palisades  Inst  for  Rsch  Services,  Inc. 
201  Varick  Street 
New  York,  NY  10014 
Attn:  Secretary 

Physics  International  Co. 

2700  Merced  Street 
San  Leandro,  C.A  94577 
Attn:  J.  Shea 

Power  Conversion  Technology,  lnc. 
11588  Sorrento  Valley  Road 
San  Diego,  CA  92121 
Attr. ; V.  Fargo 

R &£  D Associates 
P.  O.  Box  9695 
Marina  del  Rey,  CA  291 
Attn:  W.  Karzas 

C.  MacDonald 

R D -Associates 
1401  Wilson  Blvd.,  Suite  500 
Arlington,  VA  22209 
Attn:  J.  Bombard! 

Raytheon  Co . 

Hartwell  Road 
Bedford.  MA  01730 
Attn:  G.  Joshi 

Raytheon  Co. 

528  Boston  Post  Road 
Sudbury,  MA  01776 
Attn:  H.  Flescher 
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RCA  Corp.  - Government  Systen^s 
Division  - Astro  Electronics 
P.  O.  Box  300.  Locust  Corner 
East  Windsor  TouTiship 
Princeton.  NJ  08540 
Attn:  G.  Brucker 

RCA  Corp. 

David  Sarnolf  Research  Center 
P.  O.  Box  432 
Princeton,  NJ  08540 
Attn:  D.  O'Connor 

C.  Hughes 

Rensselaer  Polytechnic  Institute 
P . O.  Box  9fa5 
Troy,  .NY  12131 
Attn:  R.  Ryan 

R.  Gutmann 

Research  Triangle  Institute 
.\ttn;  Security  Ofiicer 
P.  O.  Bex  ’.2194 

Research  Triangle  Park,  NC  27709 
.A.ttn:  Eng  [>iv,  MavTant  Simons , Jr. 

Rockwell  International  Corp, 

P.  X-. . Box  5105 

-Anaheim,  C.A.  92803 

-Attn:  James  E,  Bell,  MS  HAIO 

K . Hull 
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